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PREFACE 

This reporting year we have consol i dated our four separate reports 
into one with three sections — forest inventory, forest stress, and 
standardization and cal ibration studies. A series of papers is found 
in each section. We believe this new format improves the report for 
the reader by directing him to his particular areas of interest. 

Short section summaries can be found following the listing of 

acknowledgements And cooperators . 
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SECTION SUMMARIES 


SECTION I - FOREST INVENTORY 

Microscale Photo Interpretation of Forest and Nonforest Land Classes 

A test was made to determine the accuracy of interpretation for 
thirteen forest and nonforest land classes on 1:420,000 scale high 
altitude infrared color photographs. The photographs covered eight 41- 
square kilometer (l6-square mile) study areas near Atlanta, Georgia, 
and were taken during a period of 10 months to represent three seasons -- 
early summer, late summer, and late winter. Three interpreters examined 
the photographs that were magnified 13 times using a projector-viewing 
device. On the average, the three interpreters classified 97 percent of 
all forest points correctly on winter photography. Seventy-one percent 
of the combined forest types were correct on winter photography -- 80 ^ 
of the pine, (>0% of the bottomland hardwood, and 70% of the upland hard- 
wood. Only 33% of the pine and hardwood type was correctly classified. 
This type probably should not be used at this scale. Noriforest classi- 
fication accuracy was poor but, most important to foresters, only (>% of 
all nonforest points were mi sclass if led as forest. Winter photography 
was best for minimizing these errors. Comparisons made with the results 
of two independent studies indicate that microscale infrared color photog 
raphy can be useful for forest area estimates and as a basis for strati- 
fying, -forest area. in- mul t istage forest - inventor-ies . - . - _ 

Mul tiseasonal Film Densities for Automated 
Forest and Nonforest Land Classification 

. _ M_jitiseasona.l Ektachr.ome infrared -color.-photographs taken f rom -high - 





altitude and from space were scanned using an automatic recording micro- 
densitometer. The red, green, blue, and neutral film densities were 
related to ground truth to develop microscale photo density signatures 
for thirteen forest and nonforest land-use classes. Mean densitites 
were plotted and examined for three seasons -- early summer, late summer, 
and winter. Discrimination between classes is shown to be best on 
winter photography using red density. However, standard deviations of 
the class means overlap, reducing the possibilities of discrimination 
based upon density alone. Two sources of unwanted density variation 
are examined -- film exposure variations between flights and film expo- 
sure variations caused by both d i f ferences in sun angle (season) and 
differences in film spectral sensitivity to varying levels of light. 

These sources of variation are demonstrated using mean red density for 
pine forest and pasture. Further research is recommended to isolate and 

make adjustments for this undesirable variation. 

\ 

Classification of Land Use by Automated Procedures 

Automated pattern recognition methods to identify land use were 
used on microdensitometer scans of Ektachrome i nfrared photographs for 
each of three high altitude aircraft missions and the Apollo 9 photo- 
graphic mission. Feature variables used were combinations of densities 
from red, green, blue, and clear filter scans. The results from two 
classification procedures -- linear discriminant analyses and empirical 
distribution analyses -- were compared for accuracy of classification. 
Both methods were considerably better than would be. expected by random 
assignment. In general, the empirical analyses were better than the 


2 



linear discriminant analyses. The. results from Apollo 9 photography 
were as good as those from the high altitude aircraft photography. 

SECTION II - FOREST STRESS 
Trend and Spread of Bark Beetle 
Infestations from I968 through 1971 

The spread, growth, and change of bark beetle infestations are 
dynamic features operating wl thl n the forest biome. We can follow these 
changes in the Black Hills of South Dakota on a yearly basis by counting 
tree mortality on medium scales of color aerial photography plus some 
ground checking. In four years, infestations have vacillated on a small 
area (8 square kilometers) from 97 in 1968, to 211 in 1969, to 96 in 
1970, to 2^3 in 1971 ; over 7,000 trees have been killed within these 
infestation centers. A map is presented showing the clustering and 
periodic effects of an active uncontrolled epidemic over the four-year 
period. Regressions have been computed and drawn to show the relation- 
ship of photo tree counts to ground tree counts; all of the r^ values 
were above 0.95- 

Aerial Photos Determine Optimum Levels 
of Stand Density to Reduce Bark Beetle Infestations 

Bark beetle epidemics usually get started in pine stands that are 
undergoing stress of some kind -- lack of moisture, crowding, or disease. 
Many of the pine stands in the Black Hills are too dense and need thin- 
ning. Markets for selling the wood are just now emerging. This paper 
describes how aerial photos can be used to define timber stands most 
. _l.i.kely to. succumb to beet le, attack. The-method wi ll di rect the-forest-- 
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manager to those stands needing thinning first. A crown closure-crown 
diameter comparator was used on the aerial photos to relate the Images 
to stand density conditions on the ground. On the ground, stand density 
was measured for infested stands (22), noninfested stands (16), and 
thinned stands (14) using basal area and number of trees per acre as 
criteria of stand density. We found that no recent mortality occurred 
in thinned stands and that tree diameters were significantly larger in 
infested stands. Crown closure on the photos is positively correlated 
with basal area while crown diameter is negatively correlated. 

The Use of Airborne Spectrometers 
and Multispectral Scanners for Previsual Detection 

of Ponderosa Pine Trees under Stress from Insects and Diseases 

This report outlines the work done over the past seven years to 
detect stress in conifers -- both, at Michigan in, 1964 and in the Black 
Hills from 1966 to the present. Many improvements are noted in the, air- 
borne scanning instruments, but greater advances were made in digital 
preprocessing and analog processing techniques. Despite spectral and 
radiant existence differences measurable on the ground, the accuracy of 
airborne detection is not adequate with scanners available prior to 1971. 
Indications are that a single-aperture system, such as the Bendix 24- 
channel mul t i spectral scanner or the new Michigan M-7 system, may provide 
better accuracy levels. Both the Earth Resources Program and ERTS pro- 
gram will share the same biophysical site where three 80 -foot towers 
have been erected in anticipation of a mul tispectral flight in May 1972 
and ERTS data every 18 days. 




Establishment of Earth Resources Technology Satellite Test 
Site in the Black Hills and Identification of Natural Resolution Targets 

This report describes the location of an ERTS forest and range test 
site in the Black Hills of South Dakota. Natural resolution targets 
that will be located by mercator grid coordinates include 50 replicated 
samples of each of the following: (l) healthy forest, (2) insect 

Infestations of various sizes, (3) rangeland (both natural and planted), 
(A) soil-rock outcrops, and (5) water bodies of various sizes. One 
newly infested pine stand of 5^ trees was selected for installation of 
field instruments to interface with a data collection platform. The 
location of this site is shown on a test site map. 

Detection of Root Disease Impacts on Forest Stands by 
Sequential Orbital and Suborbital Mul t i spectral Photography 

The greatest disease impact to valuable Douglas-f i r-hemlock stands 
of the Pacific Northwest is caused by Pori a weiri i root rot. To study 
this impact, three test sites were established in forested areas of 
Oregon and Washington -- two in the high Cascades and one in the Coast 
Range. Each test site is 9 square miles and represents a wide variety 
of terrain, vegetative cover, and forest and disease conditions. Color, 
color infrared, and black-and-white photographs at scales ranging from 
1:4,000 to 1:31,680 were taken of the test sites. Interpretation of the 
photographs indicates that Poria disease indicators are more accurately 
discerned on the two test sites in the high Cascades than on the test 
site in the Coast Range of Oregon. Preliminary analysis of the photo 
interpretation results indicates a range in-accuracy from to 33% on - 
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the 1:15,840 scale photographs. Accuracies of 87^ and 73% were found 
using 1:31,680 scale photographs on the two Cascade test sjtes. It 
appears that imagery taken at much higher altitudes will also detect 
the presence of the root-rot disease signatures. 

SECTION III - STANDARDIZATION AND CALIBRATION STUDIES 
Calibration of Color Aerial Photography 

In this report, the second by this research unit dealing with sensi- 
tometric calibration of aerial films, we have gone into detail concern- 
ing the spectral balance of the sensitometer or light source and the 
matching of wavelength distributions to daylight sources out to 920 . 
nanometers. A 1 i terature survey shows notable variations in data for 
different combinations of direct sunlight and scattered sky 1 ight . Several 
system spectral distributions have been calculated for possible lamp- 
filter combinations, and least-square fits to. the daylight curves were 
performed. An additional technique of electro-optical ly measuring the 
sensitometer effective exposure time and repeatability is included. 

The use of calibration techniques in testing films is further exem- 
plified in a section describing a recent flight test of new color aerial 
film emulsions. Useful film speed data were derived, and heretofore 
unobtainable details for spiral reel processing of Eastman SO-397 film 
are given. Eastman 2443 infrared color film was successfully tested 
with and without a corrective blue filter. 

Finally, a discussion of the importance of sensitometry in automatic 
interpretation of aerial film is presented. Several hypothetical exam- 
ples exhibit the probable errors introduced by working with uncalibrated 
density alone in the discriminant analysis of aerial infrared color film. 
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MICROSCALE PHOTO INTERPRETATION OF FOREST AND NONFOREST LAND CLASSES 

by 

Robert C. Aldrich and Wallace J . . Greentree, 

INTRODUCTION 

Foresters use aerial photographs to separate forest from nonforest 
land and to stratify forest land by types, volumes, and site classes, 
this is usually the first step in an extensive forest inventory such as 
the nationwide Forest Survey.^ The accuracy of forest area estimates 
and the efficiency of forest sampling designs depend upon this first 
level of information. 

Until recently, the interpretation of forest land has been done on 
prints made from panchromatic or infrared-sensitive films of scales rang- 
ing from 1:15,8A0 to 1:30,000.2 The gains from using these photographs 
could usually be shown in terms of reduced sampling errors and reduced 
inventory costs. However, two factors still limit the use of aerial 
photographs in forest inventories: (1) the increasing cost of aerial 

photography and (2) the need for more frequent flights to update resource 
information. 

^Forest Survey is a branch in the Division of Forest Economics and 
Marketing Research, Forest Service, U*S. Department of Agriculture, 
Washington, D. C. The Forest Survey was authorized by the McSweeney- 
McNary Forest Research Act of May 22, 1928. 

^Color negative films are now being used for resource photography 
oh some of the National Forests in the West. 
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If remote sensing of any kind is to show greater benefits for 
forestry, it must be up-to-date. Changes in forest area and the condi- 
tion of forest lands are occurring at an increasingly rapid rate. In 
the more highly populated areas this rate of change is becoming alarm- 
ing. To use out-of-date photographs as a basis for forest inventories 
causes a conflict between photo predictions of conditions and ground 
truth. This usually means that little gain, if any, can be shown from 
using them. 

One possible way of reducing costs and increasing gains, while 
still maintaining accuracy objectives, would be to use smal 1 -seal e or 
even microscale aerial photographs. One disadvantage of conventional 
aerial photographic scales is the large number of photographs required 
to cover an area of any size. For example, one 1:20,000 scale photograph 
( 9 - X 9“inch format) covers an effective area of approximately 5 square 
kilometers (2 square miles). A reduction in scale to 1:60,000 would 
increase the effective coverage to approximately 57 square kilometers 
(22 square miles) or over a 10-time increase in area coverage. One can 
go further and show that by using a 9" x 9“inch photograph from high 
altitude, or from space, taken at a scale of 1:400,000, we could effect- 
ively increase our coverage to approximately 2,460 square kilometers 
(950 square miles) -- over a 450-time increase. 

But what kind of information will this microscale imagery provide, 
and with what accuracy? The research covered in this report is intended 
to answer these questions and to help foresters understand how microscale 
photography might help them in their forest inventories. However, it 
should be pointed out that these results are preliminary and should be 
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accepted only on that basis. 

METHODS AND PROCEDURES 
GROUND TRUTH 

Ground truth for the Atlanta test site consists of a combination 
of several scales of aerial photography and ground observation points 
on fifteen A-mi 1 e-square study areas. In March 1970, 1:32,000 Ektachrome 
infrared color film (84^3)^’^ and 1:12,000 Ektachrome (MS) color film 
(2^»A8) , developed to a negative, were taken over each of the study areas. 
These photographs, plus 1:2,000 and 1:12,000 color taken in April 1969, 
were used to delineate thirteen land-use and forest classes. Next, 
over two hundred random points representing all forest and agricultural 
nonforest classes were selected for ground examination. Among these 
points were fifty 0.2A-hectare (O. 6-acre) forest plots taken to establish 
forest volume, forest type, and site classes. Another eighty forest, 
points were established to record forest types ,, spec i es composition, 
crown density, and stand-size classifications. Eighty nonforest agri- 
cultural points were located on the ground to record agricultural use 
at the time of each high altitude overflight. The details of data col- 
lection were reported by Langley, et al (^) and Aldrich, et al (^) . Final 
ground truth status maps were made by combining photo interpretation 
with observations made on the ground. 

' ‘^Tfcide names and ' commefcTal ehterpr i ses or pfodu“cts“ aTe ment ioned 
solely for necessary information. No endorsement by the U.S. Department 
of Agriculture is implied. 

^This film will be referred to during the remainder of this report 
as IR color. 
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PHOTOGRAPHY 


At the outset of this study, plans, were laid to obtain high alti- 
tude photography for the test site at four seasons of the year. Photog- 
raphy was planned for early June to represent early summer phenol ogi cal 
development when solor radiation is at its peak. Other flights were 
planned for late summer, fall, and winter to represent different periods 
of vegetation development. Unfortunately, this planned program could 

not be carried out entirely. Aircraft mechani cap problems and unfavor- 

’ / 

albe weather prevented a flight scheduled for late October. Thus, only 

three seasons of photography were obtained. The photographic sensor 

data for these three flights are shown In Table 1. Photographic coverage 

/ 

for the 15 study areas is shown in Table 2. 

As a result of poor photographic coverage for primary study blocks 
1 and 3, we were forced to modify our original plans. Blocks 1 and 3 
were dropped as primary study areas, and blocks 6, 8, 9, 10 and 14 were 
added (Figure 1). Although this meant a great deal of additional inter- 
pretation and data handling, it did' result in a large increase in informa- 
tion for our test. 

The interpretation test reported here was made using Hasselblad 
1:420,000 IR color photography. One major problem for interpreters was 
caused by the difference in film quality and differences in atmospheric 
conditions at the time of photography. These caused considerable varia- 
tion in quality of the film images between missions. For instance, the 
June (Mission 131) film emulsion had a weak cyan layer. This resulted 
in an overall hazy blue appearance and what appears to be a poor infrared 
response (Figure 2). The September (Mission l4l) film was quite good. 
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Table 1. Photographic sensor data by season (RB-57 mission). 



Early Summer 

Late Slimmer 

Winter 


(Mission 131) 

(Mission l4l) 

(Mission 158) 


June 

1970 

September 1970 

March 1971 

Camera and 


- 





focal length 

Film 

Filter 

Film 

Filter 

Film 

Filter 

Zeiss, 12" 

Ekta IR 
(SO-llT) 

W15 

Ekta IR 
(21+1+3) 

WI5 ' 

■ Ekta IR ■ 

(2I+U3) 

WI5 

Wild, 6" 

Ekta 
(SO- 397) 

2E 

Ekta 
(SO- 397) 

2E- • 

Ekta 

(SO- 397) ' 

2E ■ 

Wild, 6" , 

Ekta 
(SO- 397) 

450 ym 

Ekta IR 
(21+1+3) 

WI5 

Ekta IR 

(21+1+3) 

WI5 

Hass 5 UO mm 

Ekta IR. 
(SO-117) 

W15 

Ekta IR 

(SO-117) 

• WI5 

Ekta IR 

(SO-117) 

WI5 

Hass , UO nun. 

Pan 

(21+02) . 

25A 

Pan 

(21+02) 

25A 

Pan ■ 

(2I+02) 

25A 

Hass, ^0 nun 

Pan 

(21+02) 

58 

Pan 

(2I+02) 

58 

Pan 

(2I+02) 

• • 58- 

Hass , UO mm 

Pan 

(21+02) 

W12 

Ekta IR 

(SO-117) 

WI5+ 

CC30B 

Ekta IR 

(SO-117) 

W1.5+ 

CC30B 

Hass, ^0 mm 

IR 

(21+21+) 

89B 

IR 

(2I+2I+) 

89B 

IR 

(2421+) 

89B 

Hass , 1*0 mm 

IR 

(21+21+) 

W12 

Ekta IR 

(SO-117) 

WI5+ 

CC30B 

Ekta IR 

(SO-117) 

WI5 


In this table, and elsewhere in our report, "Hass" is used as the abbreviation 
for "Hasse 1 blad . " 
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Table 2. Photographic coverage for five primary study blocks and ten 
secondary study blocks by season (RB-57 mission). 


Mission 

Season 

Date 

Primary Blocks 

12 3 4 5 

Secondary Blocks 
6789 10 11 12 13 14 

131 

Early 

Summer 

June 8, 
1970 

X 

X 

1 

X 

X 

X 

X 

X 

X 

i 

X 

i 

X 

X 

X 

X 

ll+l 

Late 

Summer 

Sep. lU, 
1970 . 

1 

1 

1 

X 

X 

X 

1 

X 

X 

X 

X 

X 

1 

X 

15.8 

Winter 

March 5 » 
1971 

1 

X 

X 

X 

X 

1 

1 

X 

L_ 

X 

X 

1 

X 

1 

X 


15 

X 

















































-Figure 1. The- Atlanta test site includes all or portions of- 27 counties in 
Alabama and Georgia. Eight intensive study areas used in this report are 
shown with heavy hatch marks and solid boundaries; Two additional intensive 
study areas (hatch marked) were dropped because of insufficient high-alti tude 
photographic coverage. The remaining study blocks (secondary) will be used to 
test Interpretation models in another phase of the study. 
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Figure 2. These three 1:420,000 scale IR 
block 14 represent the quality of imagery 
June 8, 1970 , (B) September 14, 1970, (C) 
and two nonforest classes are pointed out 
hardwood, (3) pasture, (4) abandoned, and 


color photographs for study 
used in this study: (A) 

March 5, 1971. Three forest 
in (B) : (1) pine, (2) upland 
(5) bottomland hardwood. 
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and the Infrared response was enhanced by the addition of a CC+30B color 
correction filter. This same enhancement filter was used with the 
March (Mission 158) film, but there was a much more apparent vignetting 
effect causing the centers of the pictures to be overexposed. This 
resulted in some problems in resolving low-contrast details during inter- 
pretation. 

INTERPRETATION 

Three interpreters examined the 1:420,000 IR color photographs for 
all primary study areas covered by each mission. One of these inter- 
preters was a student in biological science with no previous experience 
in photo interpretation. The other interpreters were experienced and 
had worked in both the ground and photo interpretation phases of the 
study for two years. This experience does not bias the study results 
because of the interpretation techniques that we used. 

Each interpreter was given a period of indoctrination and training 
prior to beginning the test. Two study blocks (blocks 1 and 3) not used 
in the test were used for training purposes. The interpreters were 
allowed to examine all available photography including the 1:32,000 and 
1:420,000 IR color, and also were provided with ground truth to make cor- 
relations for interpretation purposes. 

Because images on 1:420,000 scale imagery are much too small to 
interpret by normal methods, they were enlarged approximately 13 times 
using a projection-viewer designed for this purpose (Figure 3). The 
device was constructed using a Bell and Howell 35 mm slide projector, an 
adjustable mirror, and a Polacoat Lenscreen viewing surface. Enlargements 
of from 8 to 22 magnifications and adjustments for tip and tilt are 
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Figure 3- This projection-viewer was used to enlarge 1:420,000 scale 
photographs to coincide with 1:32,000 ground truth strip maps; (A) Bell 
and Howell 35 mm projector, (B) adjustable mirror, (C) adjustable view- 
ing surface, (D) focusing adjustment, and (E) pulley for adjusting mirror 
distance for scaling purposes. 
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possible using the instrument as i t is presently designed. 

Each transparency was mounted as a x 4- inch lantern slide 

for insertion in the projection-viewer. The position of the slide could 
be adjusted to project the portion of the area to be interpreted. 

Interpreters, working independently, outlined and classified the 
13 land-use and forest types along each of 18 sample strips. Inter- 
pretation always began with the June 1970 (Mission 131) imagery arid pro- 
gressed to the September 1970 (Mission l4l) and March 1971 (Mission 158) 
imagery in that order. The appropriate slide was inserted in the pro- 
jector and enlarged so that strip beginning and end points, scaled from 
the ground truth strip maps, coincided with the end points on the pro- 
jected images. Interpreters were allowed to use magnifying glasses and 
filters (Wratten 15 and 25) whenever desired to enhance separations 
between some classes. The land classes were carefully mapped for each 
strip on acetate overlay material. 

ANALYSIS 

Random points selected from the ground truth records were used to 
check the photo interpretation. The number of points in each land-use 
class by season is shown in Table 3. 

Templates were made up for each strip to show the center of each 
land-use class selected to check photo interpretation. These templates 
were drawn to scale from the ground truth strip maps. Next, the templates 
were laid on top of the overlays made by each interpreter and the land- 
use or forest classification at each point noted and recorded. 

Because of distortions inherent in small-scale photographic imagery, 
and because of limitations in projection systems, we made certain 
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Table 3. Number of random check points in each, land-use class by season of 
photography. 



Leind Use^ 


Season 

1 

2 

3 

k 

5 

6 

, 7 

8 

9 

10 

1 1 

11 

12 

13 

Total 

Early 
Summer 
(June ) 

69 

37 

4o 

54 

5 

23 

30 

10 

.35 

1 

22 

13 

1 

343 

Late 
Summer 
( Sept ) 

59 

33 

l6 

45 

7 

8 

20 

6 

22 

0 

18 

7 

0 

24i 

Winter 

(March) 

6o 

28 

36 

4t 

i 

15 

23 

28 

9 

29 

3 

21 

11 

1 

311 


1 - P i ne 

2 - Pine-hardwood 

3 - Bottomland hardwood 

4 - Upland hardwood 

5 - Crop 

6 - P 1 owe d field 

7 - Pasture 

8 - Idle 

9 - Abandoned 

10 - Orchard 

1 1 - Urban 

12 - Turbid water 

13 “ Clear water 
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concessions to the interpreter. If the classification, as shown by the 
interpreter was on or within 1 millimeter of the correct ground classi- 
fication (on the enlarged photo), the interpreter was considered correct 

Once all data had been recorded they were punched on IBM cards for 
computer analysis. Computer tabluations showing the frequency of inter- 
pretation by ground classifications were used to analyze the accuracy 
of interpretation and to point out sources of mi sclass if i cat ion . Becaus 
the number of observations in some classes was too small for analysis, 
we have combined them with the land use most closely associated. Thus, 
class 5 was combined with class 6, class 8 with class 9, and class 13 
wi th class 12. 

RESULTS 

While considering the results shown below, the reader should keep 
in mind that the photo interpretation was done on 1:420,000 scale photog 
raphy. At first glance the results look rather poor. However, a com- 
parison between these results and the results of studies made in the 
past using conventional 1:20,000 scale photography shows some good cor- 
relations. This comparison is made in the CONCLUSION AND DISCUSSION 
section of this report. 

LAND-USE CLASSIFICATION 

Forest land was separated from nonforest land classes with better 
than SS% accuracy (Figure 4). As might have been expected the agri- 
cultural uses were least accurately identified under this land-use class 
ification scheme. The category "active agriculture" (crops and plowed 
fields) was classified correctly most often on the June photography, but 
even then the interpreters were correct only 51^ of the time, on the 
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Figure 4, The average accuracy is shown for three interpreters classifying 
land use on 1:420,000 IR color photographs taken during three seasons. 

Note that only in crops and plowed fields Is there a significant difference 
that might be attributed to season, June, or early summer photography, appears 
to give the greatest accuracy. 
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average. On March photography the accuracy dropped to a low of 19^- 
Pasture land is interpreted correctly in 90^ of all chances on March 
photography. On June and September photography the accuracy drops to 
Th% and 1S%, respectively. Regardless of the season, idle and abandoned 
land is mapped correctly in approximately 30^ of the cases. Orchards 
(pecan and peach) are correct in fewer than 8% of the total number of 
chances . 

Urban land, including improved roads, highways, power lines, pipe 
lines, and land areas in and around communities and cities not qualify- 
ing as forest or agriculture, is classified correctly over 97% of the 
time on June photography. The accuracy in September was In March 

the accuracy fell to 87% because some improved roads and highways could 
not be seen against highly reflective backgrounds. 

Water, including small farm ponds, can be detected 91% of the time 
in September (Figure 4). Less accuracy occurs on June and March photog- 
raphy primarily because of poorer definition due to overexposure of center 
portions of the photographs. 

There were very few forest points mi sclass i f ied as nonforest (Figure 
5). The season of photography had little effect on the results, and 
there is little indication that the classification errors can be attrib- 
uted to anything more than random interpreter error. By taking more 
care and by improving our training materials we could expect to reduce 
these errors to almost zero. 

Nonforest points were called forest quite frequently. This is a 
particularly serious problem, in the idle and abandoned and orchard land 
classes. Idle and abandoned land appears a dark gray tone and moderately 
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Figure 5- This chart indicates the percentage (average of three inter- 
preters) of all nonforest points that were mi sclass i f ied by land use and 
season of photography. Note that active cropland and plowed fields and 
the idle and abandoned land are most often mi sclass i f ied as pasture land 
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rough in texture, very much like low-stocked upland hardwood. This 
similarity is caused by clumps of weeds, blackberries, and noncommercial 
tree species scattered around with accumulated dead plant debris. It 
is easy to overestimate the stocking of forest trees and misclassify 
the land as a result. Another problem is orchards. Pecan and peach 
orchards look very much like upland hardwoods or abandoned agriculture. 
For this reason, most orchards were called either forest or idle and 
abandoned land. Orchards are extremely limited in the test areas, and 
as a result we had only a small number of samples. It is possible that 
with more extensive areas in this class, we could have shown better 
results. 

There were very few miscl ass i f i cat ion errors in the urban and water 
classes. The most frequent cause of errors in urban classification was 
small wooded areas (green belts) or garden crops within urban centers 
(towns). Another cause of mi scl ass i f i cat ion was poor contrast between 
roads and surrounding agricultural land caused by overexposure, particu- 
larly on the March photography. Water was mi scl ass i f i ed in almost every 
instance because of the poor contrast which (turbid) water has with sur- 
rounding classes, particularly in overexposed portions of the photographs 

FOREST CLASSIFICATION 

The accuracy of forest classification varies considerably by type 
and by season of photography (Table 4). Photo interpreters are shown to 
be relatively consistent with one possible exception. Interpreter A was 
inexperienced and apparently had insufficient training in relating the 
differences in infrared response on March photography to vegetation types 
Despite this, we conclude that photo interpretation on March photography 
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Table 4, Accuracy of forest type classification by photo interpreter and 
by season. 


! 

1 



Forest 

Type^ ’ ^ 



Photo 

. 1’ 

2 

3 . 

4 

Season^ 

Interpreter 

(P) 

(PH) 

(BH) 

(UH) 



— — 

- - percent 

correct - - 

— — 


A 

71 

0 

55 

6l 

June 

B 

70 

11 

25 

72 


C 

67 

0 

48 

71 


Me an 

69 

4 

^,3 

68 





A 

79 

12 

6 

74 

September 

B 

83 

15 

13 

67 


C 

1 

8l 

' 33 

6 

' 71' 


Mean 

81 

20 

8 

70 



1 - 


A 

67 

7 - 

53 

55 

March 

B 

80 

28 

64 

79 


C 

82 

39 

64 

• 64- 


Mean 

76 

25 

6o 

66 


^Seasons tested: 

June - early summer 
September - late summer 
March - late winter 


2 

Forest type: 

Pine (P) - 1 
Pine-hardwood (PH) - 2 
Bottomland hardwood (BH) - 3 
Upland hardwood (UH) - 4 


^The number of observations by forest type is shown in Table 3- 
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resulted in the best forest type classification. Interpreters B and C 
were able to classify pine type correctly 8l| of the time. However, 
these two interpreters could classify pine-hardwood type correctly only 
33% of the time; this low level of accuracy is not unusual even on 
conventional 1:15,840 and 1:20,000 scale photography. Part of the dif- 
ficulty is in defining the percentage of pine in the stand. ^ We found 
that with the better resolution of 1:12,000 and larger scale color and 
infrared color photography this classification becomes easier and more 
accurate to identify. 

The greatest advantage of high altitude photography taken in March 
is found in the greater accuracy of hardwood type classifications. On 
the average, 60^ of the stands called bottomland hardwood and 66% of the 
stands called upland hardwoods were classified correctly (Figure 6). In 
contrast to this, only k3% of the bottomland hardwood was correctly class- 
ified on June photography, and only 8^ on the September photography. 

The accuracy of upland hardwood interpretation does not appear to vary 
significantly with the three seasons tested. 

There were very few forest points mi sclass i f led as nonforest land 
(Figure 7). Regardless of the season, no fewer than 96^ of forest points 
were correctly classified as forest. Thus, errors in forest type classi- 
fication were largely errors in judgment of stand stocking. For instance, 
on the average, interpreters classified 19^ of the pine stands as upland 

^Pine- hardwood: 25“50^ of the dominant stand is in pine; the 

remainder is composed of hardwood (deciduous) tree species. 
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Figure 6. This chart shows the average accuracy for three interpreters 
classifying forest types on 1:420,000 IR color photographs taken during 
three seasons. Note that forest types are generally most accurately 
identified on March photography. 
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hardwood on June photography. On September photography the error was 

and on March photography the error was 14%. These are relatively 
consistent errors and should be considered the most serious errors in 
forest misclass i f icat i on. These errors may be reduced in the future 
by improving photo resolution and developing better keys and definitions 
of forest types. 

The greatest errors in pine-hardwood type classifications were in 
calling pine-hardwood mixtures upland hardwood (Fig. 7). Sixty-two 
percent of the pine-hardwood stands were called upland hardwood on 
June pho.tography , 53% on September photography and 43% on March photo- 
graphy. It is apparent that as the percentage of pine is reduced in 
mixed stands, the chance for calling the stand hardwood is increased. 
Also, when deciduous trees (hardwoods) are leafless, there is a greater 
chance of seeing pine in mixed stands and a greater chance of calling 
the stands correctly. 

Bottomland hardwood type was often called upland hardwood type 
(Figure 7). On the other hand, only occasionally (average of 11^) 
was an upland hardwood stand called bottomland hardwood. Mi sclass i fi ca- 
tion of bottomland hardwood as upland hardwood was greatest in September 
and least in March — 77% as compared to only 24%, Therefore, when we 
want to distinguish between upland and bottomland types, winter photo- 
graphy is most useful. 

CONCLUSION AND DISCUSSION 

It is always difficult to evaluate the results of a study such as 
this without some basis for comparison. For instance, the results shown 
here were arrived at using 1:420,000 scale infrared color. Do we know 
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Figure 7- The percentage (average for three interpreters) of all forest points 
that were mi sclass i f ied by type and season of photography. Note that pine-hard 
wood and bottomland hardwood types are most frequently misclass if ied. In both 
cases they are usually called upland hardwood by interpreters. 
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how well Interpreters might do on conventional 1;20,000 scale panchromatic 
photographs? 

A search of the literature shows only two studies that can be used 
for comparison. In the first (_][), interpreters used basically the same 
set of criteria to classify land use in southwest Georgia on 1:20,000 
panchromatic photographs. In this example 95.5% of all forest classi- 
fications were correct. This compares favorably with the 96% or better 
accuracy shown in this report. The second study taken from the litera- 
ture shows the expected accuracy for forest type classification on 
1:20,000 scale photography Q) • This test made by the TVA in 1952 showed 
that interpreters could classify three types (pine, mixed, and hardwood) 
correctly 7^% of the time. When the forest types used in the present 
study are combined in a similar way, we can show that 71% of the type 
classifications were correct if made, on March photography. 

These two comparisons can hardly be considered conclusive evidence 
of the benefits of 1:420,000 scale infrared color photography for forest 
and nonforest classification. However, they should lend some credence 
to the statements below. 

1. Microscale IR color photography (1:420,000) can be interpreted 
within reasonable limits of error to estimate forest area. 

2. Forest interpretation is best on winter photography, and provides 
97% or better accuracy. The greatest source of error is in calling 
abandoned agricultural land forest land, but this error is minimized on 
winter photography. 

3. Broad forest types can be classified on microscale photography. 
For instance, pine type is correctly identified more than 80% of the time 
on winter photography and at least 70% of the time regardless of season. 
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Pine-hardwood cannot be correctly identified better than 25% of the time 
(winter photography) and probably should not be attempted on microscale 
photography in the future. Bottomland hardwood is correctly classified 
60% of the time on winter photography but cannot be separated from 
upland hardwoods at other seasons of the year. Upland hardwood is usually 
classified correctly 70 % of the time regardless of season. 

4. Active agricultural land is classified most accurately on early 
summer photography; in this test 51% of active cropland and plowed fields 
were correctly classified. The most frequent error was interpreting 
cropland as pasture. Fortunately, only a very small percentage of the 
cropland category was misc lass i f ied forest land, and on winter photography 
this error was at a minimum. 

5. Only 6% of all nonforest observations (including urban and 
water) were misclass if ied as forest. Winter time was the best season 
for minimizing these errors. 

During the coming year we will continue photo interpretation tests 
using high altitude aerial photography. These tests will include rigid 
statistical tests of the data presented here. For instance, a weighting 
scheme has been devised to assign to misclass i f ied plots -- the magnitude 
of the weight depending on the seriousness of the error. An analysis of 
variance will be made for these weighted data to determine the levels of 
significance for interpretation errors. 

Mu 1 1 i seasona 1 color and IR color photography taken at a 1:120,000 
scale will be interpreted and analyzed in a manner similar to that used in 
this report. We will also concentrate on developing a photo interpretation 
key for microscale photo interpretation of forest and nonforest classifi- 
cations. 


30 



LITERATURE CITED 


1. Aldrich, R. C. 1953. Accuracy of land-use classification and area 
estimates using aerial photographs. Jour. Forestry. 51(1):12-15. 

2. Aldrich, R. C., W. J. Greentree, R, C. Heller, and N. X. Norick. 

1970 . The use of space and high altitude aerial photography to 
classify forest land and to detect forest disturbances. Annual 
progress Report for Earth Resources Survey Program, OSSA/ 

NASA, by the Pacific Southwest Forest and Range Experiment Station, 

36 p., illus. 

3 . Bsteson, A. R. 1952 . A test of aerial photography for forest inven- 
tory purposes in the Tennessee Valley. Report No, 205 - 52 . Forestry 
Investigations Branch, Tennessee Valley Authority, 3^ p. 

4. Langley, P. G., R. C. Aldrich, and R. C. Heller. 1969. Multi-stage 
sampling of forest resources by using space photography -- an Apollo 9 
case study. Second Annual Earth Resources Aircraft Program Status 

I Review. Vol . 2. NASA Manned Spacecraft Center, Houston, Texas. 

28 pp. illus. 


31 



MULTISEASONAL FILM DENSITIES FOR AUTOMATED 


FOREST AND NONFOREST LAND CLASSIFICATION 

by 

Wallace J. Greentree and Robert C. Aldrich 
INTRODUCTION 

This is the second of two reports on research started in 1969 to 
separate 13 forest and nonforest land-use classes using film density. 

The first report in October 1970 (1) was based upon fragmentary data for 
one study area at one season of the year. The continued research reported 
here is based upon data for eight study areas using photography for three 
seasons . 

The objectives of this research are: (I) to develop microscale photo 
density signatures for forest and nonforest land-use classes to use in 
multivariate analysis techniques, (2) to examine multivariate analysis 
techniques to find the most effective techniques for classifying forest 
and nonforest land use, (3) to compare film densities by season and iso- 
late sources of variation not caused by land use, and (4) to examine 
alternatives for removing undesired sources of variation to improve forest 
and nonforest discrimination. Research involving multivariate analysis 
techniques is reported by Norick and Wilkes in another section of this 
report. However, the association between the basic data used in these 
two reports should always be kept in mind, 

METHODS 

To meet our study objectives we used microdensitometry methods to 
evaluate multiband photographs. This involved microimage density measure- 
ments using a Photometric Data Systems microdensitometer coupled with a 
Data Acquisition Systems digitizer, and computer processing with a Univac 
1108 computer. 
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PHOTOGRAPHY 


Microscale (1:420,000) transparencies taken by NASA's RB-57 high- 
altitude aircraft were used for each of 3 seasons. The photography was 
taken with six Hasselblad 70 mm cameras equipped with 40 mm focal length 
lenses. Film and filter combinations used in this study are shown by 
mission date in Table 5. 

The Wratten 15 and 30B color correction filter combination was used 
on Missions l4l and 158 because we hoped it would be better for cutting 
haze and would improve the infrared response of infrared color film at 
high altitude (4). Our hopes were justified by the results which showed 
a much improved image sharpness and infrared response. We did not use 
this film-filter combination on Mission 131, and as a result the film 
looked hazy and had an overall blue cast. 

Master copies of 1:2,400,000 70 mm infrared color frames 3740 and 
3741 from the Apollo 9 SO-65 photography experiment were also used. 

MICRODENSITOMETRY 

Each of the 70 mm transparencies was run on the microdensitometer 
using block scanning techniques (_1^, 2 ). The density readings along our 
sample lines were then determined by computer programming. 

Each sample line on the infrared color photography was scanned in 
turn with a clear, red, green, and blue microdensitometer filter.^ The 
panchromatic photographs and the black-and-white infrared photographs were 
scanned with only a clear filter. 

ANALYSIS 

The mean densities and their standard deviations were computed for 
each of the 13 land-use classes along each random sample line within the 
sample blocks. These statistics were also computed for blocks and the sum 

^92 - red; 93 - green; 94 - blue; clear - no filter. 
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of all blocks for each RB-57 mission. Density signatures for the 
panchromatic and black-and-white infrared photographs, processed with 
a clear filter, were obtained in the same manner. The combined block 
mean density values were then plotted for visual analysis (FigureS ). 

From these graphs it was evident that some of the density variation 
was the result of sources other than land use. By intuitive reasoning 
we felt that one source might be photographic exposure. Thus, to demon- 
strate the effect of exposure differences between photos, we plotted the 
mean red density values for pine and pasture by blocks (Figure 9)- In 
theory, if the red, green, and blue filter are subtracted from the mean 
density values of the clear filter, which is unfiltered white light, we 
should, in effect, be normalizing for exposure differences. These 
values were plotted in Figure 9 for comparison-wi th the unnormalized data. 

Normalized means were computed for all forest and nonforest land-use 
classes and plotted for comparison in Figure 10. An examination of both 
the photographs and the data in Figures 9 and 10 suggests the possibility 
that density is affected by position of the study blocks within the indi- 
vidual photographs. To estimate the effect of this variable we made a 
first attempt by measuring the distance and the azimuth from the center 
of the photograph to the center of each block. The mean red density was 
plotted for pine and pasture by distance and by quadrant location (Figure 11) 

RESULTS 

Table 6 shows the number of density data points stored on tape for 
each of the 13 land-use classes by season. The number of data points in 
each column should be multiplied by 4 to include all points for clear, red, 
green, and blue density signatures. The total number of data points for 
both the panchromatic and black-and-white infrared density signatures 
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Figure 8. These graphs show a three-season comparison of mean optical red, green, and blue film densities 
for 13 land-use classes. The measurements were made on microscale infrared color transparencies (SO-117): 

(A) 1;420,000 scale taken in June 1970, (B) 1:420,000 scale taken in October 1970, (C) 1:420,000 scale taken 
in March 1971> and (D) 1:2,400,000 Apollo 9 taken in March 1969. Notice the magnitude of standard deviations 
plotted for the means in (C) . Also notable is the similarity between the plotted data in (C) and (O) . 







> 3 : 3 a 9 tc 

I 1 : 

'f summer 




Figure 9. Mean red density is plotted by study block for pine and pasture 
to show the effects of season on film density: (A) early summer, (B) late 

summer, and (C) winter. The upper portion of each graph shows unadjusted 
red density; the lower portion of the graph shows normalized red density 
(red density - clear density). Notice that the scatter of density points 
around the horizontal line of equal density has been reduced considerably. 
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RED DENSITY SIGNATURE | GREEN DENSITY SIGNATURE 
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with densities on panchromatic (C) and black-and-white infrared (D) . The seasonal responses are shown 
and indicate that red film density on winter infrared color photography has the most information. However, 
density signatures on panchromatic photography would appear to be able to separate forest from nonforest 
with greater likelihood of success. 






c 


Figure 11. The mean density measured on infrared color film is plotted 
for pine and pasture by season, distance from the photo nadir, and quadrant 
location: (A) early summer, (B) late summer, and (C) winter. The solid 

1 i ge (l) in each graph shows data located in the two northernmost quadrants 
(0 - 90 °, 270°“360°) . The dashed line (ll) in each graph shows the data 

located in the two southernmost quadrants (90°-l80°, l80°-270°) . Notice 
the general increase in density with distance and the higher density for 
southern quadrants on winter photography (low sun angle). 
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Table 5. Film and filter comb i na t i ojis used on the various RB-57 high- 
altitude photographic missions. 


RB-5T 

Mission 

Season 

Film 

Filter 

131 

(June 8 , 
1970) 

Early 

Summer 

Ektachrome Infrared, SO-117 

Panchromatic, 2l0l 

Bla,ck & White Infrared, 2h2k 

Wratten 15 
Wratten 25A 
Wratten 89B 

111 . 1 

(Sep. ll, 
1970) 

Late 

Summer 

Ektachrome Infrared , SO-llT 
Panchromatic , 2^0^+ 

Black Se White Infrared, 2^2^ 

Wratten 15+30B 
Wratten 25A 
Wratten 89B 

158 

( March 3, 

1971 ) 

Winter 

Ektachrome Infrared, SO-117 

Panchromatic , 2 i+ 0 U 

Black 8c White Infrared, 2h2h 

Wratten 15+30B 
Wratten 25A 

Wratten 89B 

1 













Table 6. Number of density data points by land use by season of photography; 
RB-57 and Apollo 9 infrared color. Combined study blocks,* 


TOTAL NUMBER OF DENSITY DATA POnfTS 


Land-\ise 

Class 

Early 

Summer 

Late 

Summer 

Winter - 

Winter 
Apollo 9 

1 

1844 

1527 

1554 

3133 

2 

300 

206 

229 

393 

3 

376 

187 

372 ! 

674 

4 

819 

552 

703 

997 

5 

75 

32 

77 

110 

6 

126 

36 

127 

155 

7 

866 

665 

883 

1174 

8 

52 

39 

53 

59 

9 

282 

210 

289. 

423 

10 

24 

0 

25 

33 

11 

518 

483 

531 

645 

12 

4l 

18 

35 

43 

13 

3 

0 

3 

3 

Total 

5326 

3955 

4881 

7842 


Early summer includes blocks 2, 4, 5, 6, 8, 9, 10, 14. 

Late summer includes blocks 4, 5 (less one strip), 6, 8, 9, 10, 14, 
Winter includes blocks 2, 4, 5, 8, 9, 10, 14. 

Winter (Apollo 9) includes blocks 1, 2, 3, 4, 5, 6, 8, 9, 10, 14. 
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is approximately the same as those shown in the table. 

The combined block mean density signatures for the land-use classes 
are shown in Figure 8 . It would appear from these graphs that a separation 
can be made between most forest and nonforest classes. However, the stan- 
dard deviations about the class means overlap one another, even though - 
the mean densities are different.^ This makes It difficult to separate 
the individual forest classes or nonforest classes from each other by 
using raw density alone. This result bears out the conclusions of our 
1970 report which was based upon only fragmentary data (J_) . 

Some of the overlap in standard deviations may be the result of 
unwanted variation that can be controlled. To explain this possibility, 
density signatures for pine and pasture were plotted by block, and the 
trend of the data was examined for possible sources of block variation 
(Figure 9 )* The plotted data shows the effects of exposure differences by 
season of photography. Ideally, the data would form a horizontal line of 
equal densities. Any differences in level between lines would be caused 
by variations in phenological development of the vegetation, variations 
in site, and differences in solar radiation at the three periods of the 
year. However, this is not the case, and deviations from the horizontal 
line of equal densities in the example must be attributed to one or more 
of the following: ( 1 ) variations in reflected energy received at the 

sensor, ( 2 ) differences in film exposure, ( 3 ) differences in film sensi- 
tivity within and between emulsions, ( 4 ) differences in spectral sensi- 
tivity of the film layers to differing light levels and ( 5 ) the position 
of the block within the photograph format, i.e., distance from the nadir. 

^See Figures (C) for standard deviations on winter photography. 
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Figure 10 Is a three-season display of density signatures for all 
land-use classes after normalizing for exposure differences. Normaliz- 
ing consisted of subtracting the mean clear density (white light) from 
the mean red, green, and blue densities on the infrared color. 

Normalizing reduced the difference between winter and late summer 
red density signatures but had little effect on reducing early summer 
red density. The difference in infrared response is shown to be best 
on winter photography. Late summer and early summer red densities show 

I 

fewer possibilities for discrimination even after normalizing. However, 
there are some interesting reversal patterns shown by season (Figure 10). 

For instance, crops, idle land, and urban all show red density reversals 
on early summer photography when compared with late summer photography. 

Clear water on winter photography shows a density reversal when compared 
with early summer photography. These patterns and others may be useful 
for discrimination of targets that cannot be separated on any, one season 
of photography alone. 

After the green and blue film densities were normalized, the density 
differences between land classes were reduced (compare Figure 10(B) with 
Figure 8), and the apparent seasonal differences that existed before 
normalizing are almost gone. From this there appears to be very little 
value in the blue and green densities. However, computer discriminant 
analysis techniques may find that differences between single-band densities, 
by season and between seasons, are valuable for separating some classes. 

The mean densities for combined blocks indicate that forest classes 
can be discriminated from nonforest classes on panchromatic (25A) 
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photography (Figure lO), in fact, even when standard deviations are consi- 
dered, there appears to be a greater potential for discrimination on this 
film than on infrared color. Of the two seasons for which photo coverage 
was available, winter was the best. Unfortunately, there is still very 
little possibility of distinguishing between forest types or nonforest 
classes. 

Black-and-white infrared film (89B) alone has little value for land- 
use classification unless it is used together with panchromatic (25A) and 
infrared color film. Then it would provide some classification signatures. 
Some interesting density reversals are evident between seasons. However, 
the full potential of these must be evaluated by multivariate computer 
routines. 

An important source of unwanted variation in film density is demon- 
strated in Figure II. In this figure the normalized mean red film density 
has been plotted by blocks for pine and pasture for three seasons. Densi- 
ties have been plotted by distance from the center of the photo and by 
north and south quadrants. Although the number of data points is limited 
there is some trend in the data that should not be ignored. For instance, 
quadrant location (northerly or southerly) appears to have little effect 
on the mean density in early summer. However, there is a slight upward 
trend to the data as distance from the center of the photo increases. 

This increase is not as dramatic as that shown for late summer and winter 
photography. On late summer photography the upward trend of density for 
pine is most apparent in the northern quadrants. The downward trend for 
the two southern quadrants is meaningless with only two data points, even 
though this could indicate a difference between the two quadrants. Winter 
photography shows the best separation of data in the north and south quadrants 
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and should give impetus to further tests -- particularly for this season. 
This is because winter photography shows the greatest all-around poten- 
tial for forest and nonforest land classification on microscale photography, 

CONCLUSION 

The use of film density to discriminate between forest and nonforest 
land-use classes entails many problems which must be resolved before it 
will be an operational technique. In another section of this report 
Norick and Wilkes show no better thatn 78 % accuracy in classification 
using film density in multivariate discrimination techniques. This is 
rather discouraging because human interpreters can separate forest from 
nonforest with an accuracy of 96% or better on the same microscale photo- 
graphs (see Aldrich and Greentree in the previous section of this report). 

Low accuracy using film density can be attributed to several undesir- 
able sources of variation. Two sources that may be partially controlled 
have been demonstrated in this report; these are also discussed on pages 
170-178 in the paper by Dana and Myhre. One Is the effect of exposure 
difference between and within the photographs and between photographic 
missions. This can be partially removed by subtracting clear density 
(white light) from raw density measured in the red, green, and blue bands. 
This has a normalizing effect on the data. 

The second form of unwanted variation remaining results from dif- 
fering sensitivities of the one or more emulsion layers in the film to 
variations in light intensity and quality. These variations are due to 
time of day, season, and camera optics, and occur within zones. Some 
of this zonal variation might be removed using corrections based upon 



the location of sample areas within the photographs. For instance, it 
was demonstrated that the density of the yellow and magenta layers (red 
density) increases gradually in a radial direction from the photo nadir. 
Generally speaking, film densities are also higher in the southern quad- 
rants on winter and late summer photography due to the low sun angle 
at these two seasons of the year. Thus, distance and sun angle could 
probably be programmed with density to increase the likelihood of success 
in automated photo interpretation. 

DISCUSSION 

During the coming year we will continue our study of film density 
for discriminating land-use on high-altitude aerial photography. It is 
hoped that additional imagery from RB-57 mission 191 flown on November 12, 
1971 , will be of good quality and will complete the seasonal coverage for 
the Atlanta test site. This year we will also examine 1:120,000 scale 
RC-8 photography with the intention of using both density and spatial 
frequency to increase the accuracy of classification. 

As a result of preliminary work reported here, we intend to make a 
more complete study of the effect of distance from the photo nadir and 
photo quadrant on film density signatures. It is our hope to develop 
functions for the correction of these exposure-related variations in 
film density and to incorporate these in multivariate discriminant 
analysis techniques. 
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CLASSIFICATION OF LAND USE BY AUTOMATED PROCEDURES 

by 

Nancy X. Norick and Marilyn Wilkes 
INTRODUCTION 

Pattern recognition techniques include a variety of mathematical and 
statistical methods aimed at identifying particular objects of interest -- 
and usually automatically. Each may be optimum under certain conditions. 
One of the most difficult tasks is finding what features, variables, or 
measurements are best to use for finding differences between population 
types, e.g,, forest and nonforest classes. The second is finding what 
classification techniques provide the best separation of the population 
types. Examples of features often used are: (T) statistical quantities 

such as means, variances, and other moments, (2) mathematical transforms 
of the data such as Fourier or Hadamard, and (3) direct untransformed 
measurements such as 3-color photo density measurements or digital output 
from a multispectral scanner. Due to time and cost restrictions on our 
data collection we were limited for this analysis to the use of direct 
measurements of image density on infrared color photos using red, blue, 
green, and clear filters. Until now we have concentrated on selecting 
among various classification procedures. For instance, we have explored 
this problem using multivariate linear discriminant analysis and various 
forms of an empirical distribution analysis. The next step will be the 
programming of a k-nearest neighbor procedure (k may be equal to 1, 2, 3, 
etc.). 

The linear discriminant analysis is the optimum procedure when the 
data have multivariate normal distributions. It is robust, i.e., it is 



nearly optimum for data which are not normal but come from continuous 
unimodal distributions. The empirical distribution may be better if the 
data are nonnormal in certain ways -- such as from multimodal distributions. 

A drawback is that it is necessary to have a large number of observations 

in the training set to adequately estimate the distribution. But, both 
of these methods may be totally inadequate if the data occur in strings. 

A method which may be better for this type of data Is nearest neighbor 

analysis. F i gure 12 shows examples of three data types and the appropriate 
distribution analysis. Our aim is to compare these three methods for our 
data (classification of land use). 

LINDA is a linear discriminant analysis computer program which we 
wrote to handle our data. It is flexible for our use but not for general 
use as it requires a programmer to modify the main routine for each 
problem. The empirical distribution program, is similar with respect to 
use. The number of cells and length of intervals can be changed for each 
run by a programmer. Both programs were written for the Univac 1108 in 

V 

Fortran. 

The linear discriminant method has been described by T. W. Anderson 
(^) . This analysis chooses that class i f i cat i on wh i ch has the highest 
probability of occurrence for the values of the variables. The probabilities 
are based upon the training sample estimates of the normal distribution 
parameters. The empirical distribution method likewise chooses that 
classification which has the highest probab i 1 i ty of occurrence; however, 
no assumptions such as those pertaining to normality are made about the 
probability distributions. Instead, the observed frequency distribution 
(often with some smoothing) of the training sample is used. Our empirical 
distribution analyses were done in several ways described below. 
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Figure 12. These three different data types may require di fferent methods for 
optimum discrimination between populations: (A) a data type for which linear 

discriminant analysis might be chosen, (B) a data type for which, empirical 
distribution ana lys is might be chosen, and (C) a data type for which the "nearest 
neighbor" analysis might be chosen. 
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SAMPLING THE DATA 


Four sets of photographic data for the Atlanta, Georgia, test site 
are used in this study. The aerial photography is 1:420,000 Ektachrome 
infrared color taken by NASA's high altitude RB-57 aircraft during three 
seasons -- early summer (June 1970), late summer (September 1970), and 

g 

winter (March 1971) • Master copies of the Ektachrome infrared photo- 
graphs taken from the Apollo 9 (SO-65) multiband photography experiment 
are also used, A full description of the photography, the photo and ground 
data acquisition techniques, and other analyses are given by Greentree and 
Aldrich in the previous section of this report and in the 1970 annual 
report (J_) • 

For each photo mission two samples of data were taken. One was used 
as the training set for the classification analyses and the other as the 
new set to test the performance of the procedure. For both sets of data 
the sampling method was the same. A random sample of ten points was 
selected from each of eighteen sample strips for each land-use class. 

When ten exceeded the number of observations for a land-use class on a 
sample strip a smaller random sample was taken. This method gave equal 
weights to all sample strips and made the number of observations per 
land use as uniform as possible. The effect of this sampling scheme was 
to give every land use a nearly equivalent chance of being identified. 

For a production analysis (for example, map making) one might choose 
another sampling method such as a completely random sample. This would 
have the effect of giving greater emphasis to classes with more observa- 
tions, Another sampling scheme could be devised to put greater emphasis 
on certain specified classes deemed more important by an investigator. 

^Mission 131, June 8, 1970; Mission I4l , September 12, 1970; Mission 
158, March 5, 1971 . 
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For each aircraft mission, information for both samples from micro- 
densitometer scans of 1:420,000 scale Ektachrome infrared color film was 
written on computer tapes. Information from 1:2,500,000 scale Apollo 
scans was recorded as well. For each observation sampled the land use 
and density values for scans using red, green, blue, and clear filters 
were recorded. Previously we had tested the effects of using the visual 
versus the clear filters for normalizing the exposure differences and 
found they were equivalent for our purposes. 

THE ANALYSES 

We tested several combinations of the original 13 forest and nonforest 
categories in our analyses. The simplest was all forest versus all non- 
forest. Another set of categories combined all forest types while all 
nonforest categories were kept intact -- a total of nine categories. 
Analyses were done also using the original 13 categories defined by 
Aldrich, et a 1 (J_) . 

The variables used to make the classifications were red minus clear 
densities, blue minus clear densities, and green minus clear densities. 

We had therefore a trivariate distribution. 

The linear discriminant analysis was done in the standard way. 

However, the empirical distribtuion analysis was done using several models. 
In the first, the density range of each variable was divided into three 
equal length intervals. This gave 27 possible cells for observations to 
fall into, e.g., an observation could be in interval 1 for R minus C, 1 
for B minus C, and 1 for G minus C; 1 for R minus C, 1 for B minus C, and 
2 for G minus C, etc. The number of training set observations for each 
land-use class in each cell was tabulated. That land-use class with the 
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largest number of observations in a given cell was the choice for the 
cell. Every test data observation that fell in that cell was classi- 
fied as being of that land-use class. One result of this method of 
choosing the interval length was that many observations fell in a few 
cells while most cells were empty or nearly empty. 

To overcome this the range of each variable was divided into three 
equal frequency intervals. This improved the accuracy of the classifica- 
tion, and because of this success we tried other empirical models with 
4, 5f and 6 intervals per variable. Each of these empirical models 
gave weights to land-use classes proportional to their sample size. 

We also tried two other weighting systems for Mission 158 (winter) 
data using 9 land-use categories and 4 intervals. The first gave equal 
weight to each category regardless of its sample size. The second 
weighting system gave greater emphasis to categories with fewer obser- 
vations. These models correctly classified more observations in the 
lightly sampled categories than the previously used models, but their 
overall performance, i.e., the proportion correct for all categories, 
was slightly poorer. Their use will depend upon the importance of the 
smaller categories. 

For the 2-category classification problem, in addition to making the 
weights proportional to the sample sizes, we ran some data with the weights 
of forest samples doubled. As we had anticipated, this improved forest 
class i f i cat i on . 

To verify that the classification procedures were having an effect 
we computed the number of correct classifications that would occur by 
chance. As a result we found that all methods were considerably better 
than would be expected by chance. 
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RESULTS. AND DISCUSSION 


The major results of this i n.ves t i gat i on have been tabulated in Tables 
7 through 22, Tables 7i 3, and 9 show the proportion correct over all 
categories for each photographic mission and each method. The propor- 
tions that would be expected if a completely random procedure had been 
used are also tabulated and are shown at the top of each table. All 
methods used are considerably better than random. 

For analyses that combined observations into 13 and 9 land-use 
classes, the empirical distribution method yielded considerably better 
results than LINDA. However, when only 2 classes were used, forest 
versus nonforest, the methods are roughly equivalent. 

One striking result is that classifications using the Apollo 9 
imagery are almost as good as those for the high flight aircraft missions. 
This may be due to the. averaging effect of the lower resol t ion space 
imagery. If so, better results might be obtained on high altitude 
photography if we were to average - dens i ty values for adjacent points, 
on the scan lines and on lines scanned above and below the original 
scan location. We will test this procedure as the data analysis continues. 

Tables 10-22 show numbers of observations by actual and predicted 
classes. From these data one can see which classes appear to be similar 
in dens i ty va 1 ues , 

FUTURE WORK 

The empirical distribution program is being modified to smooth the 
distribution curves. As a result of this, no observations will be classi- 
fied as '^undecided." Once this is completed the next classification method 
to be programmed will be the nearest neighbor method. 

Thus far, our analyses have been limited to sample blocks distributed 
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over the entire study area. An important question which we have not yet 
tried to answer is: Can training samples from one block be used to 
classify observations from another block? We hope to answer this 
question during the coming year. 

Another important area for investigation is the effect of the 
location of samples within the photo frame on class i f i cati on accuracy. 

We hope to learn whether location can be of value to improve the normaliz- 
ing procedure which has, until now, consisted only of subtracting the 
clear filter densities. 
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Table 7. Comparisons of methods for automated classification: Accuracy 

for 13 land-use classes by season of photography. 


METHOD 

Mission 131 
June 

Mission 1^1 
September 

Mission 158 
March 

Apollo S 
March 







Expected correct 
if random 

0.1052 

0.1196 

0.1041 

0.1066 

Linear Discriminant 

Analysis ! 

! 

1 

0.1853 

0.1686 

0.2288 

0.1998 

Empirical Distribution 
Analysis 





3 equal length 
intervals 

0.2250 

0.2532 

0. 2522 

0.2364 

3 equal frequency 
intervals 

O.2I08 

0.3133 

0.2688 

0.2565 

U equal frequency 
intervals 

0 . 23 hh 

0.3187 

0.2785 

0.2707 

5 equal frequency 
intervals 

0. 2866 

0 . 368 i^ 

0.3240 

0.2737 

6 equal frequency 
intervals 

0.2893 

0 . 359 i+ 

0.3284 

0.2879 
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Table 8. Comparison of methods for automated classification: Accuracy for 

9 land-use classes (all forest combined) by season of photography. 


1 

METHOD 

Mission 131 
June 

Mission l 4 l 
September 

Mission 158 
March 

Apollo 9 
March 


T-* 1 • 



r X upu X 0 X uii c u X X c u u 


Expected correct 

if random I 

0.2581 

0.2923 

0.2413 

0.2654 

Linear Discriminant 
Analysis 

0.3109 

1 

0.3178 

0.3264 

0.3963 

Empirical Distribution 
Analysis 





3 equal length 
intervals 

0.4781 

0.5279 

0.4764 

0.5095 

2 equal frequency 
intervals 



0.4737 


3 equal frequency 
intervals 

0 . 4679 

0.5526 

0.4982 

0.5036 

k equal frequency 
intervals 

0.4843 

0.5526 

0.4939 

0.5136 

5 equal frequency 
intervals 

0.5055 

0.5858 

0. 5289 

0.5077 

6 equal frequency 
intervals 

0.5130 

0.5740 

i 

0.5184 

0 . 5172 
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Table 9. Comparison of methods for automated classification; Accuracy for 2 
land-use classes (forest versus nonforest) by season of photography. 


METHOD 

Mission 131 
, June 

Mission 1^+1 
September 

Mission 158 
March 

Apollo S 

1 March 

1 ^ ^ 


J. • - -L 



jTiupuruxuiJ 



Expected correct 
if random 

0.5028 

1 

0.5000. 

. 0.5079 

0.5020 

Linear Discriminant 
Analysis 

0.7305 : 

0.7489 

0.7636 , , 

0.7571 

Empirical Distribution 
Analysis 





3 equal length 
intervals 

0.7023 

0.6599' 

0.7338 

0.7056 

3 equal frequency 
intervals 

■ 0.7155 ■ 

0.7672 

0.757i+ 

0.7^^1 

. 1 
A equal frequency 

intervals 

0.7168 

0.7597 

0.7758 

0 . 7601 

5 equal frequency 
intervals 

0.7339 

0.7650 

0.7697 

0.7^t88 

6 equal frequency 
intervals 

0. 7^+56 

0.7575 

0.7715 

0.7^T 

1 
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Table 10. Mission 158 - Color f I Im, f I Iter 15 + 30B. Empirica,! Distribution 
Analysis with 3 equal length intervals. 


ACrrUAL TYPE 


TYPE 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Total 

1 

458 

12 

28 . 

77 

JLj 

95 

16 ■ 

85 

17 

. 813 

2 

0 

D 

0 

0 

D 

D 

. 0 

_ _ i 

0 

0 

0 

3 

0 

0 

0 

n 

0 

0 

1 

0 

0 

0 

0 

4 

38 

26 

35 

70 

22 

33 ^ 

2 

73 

8 

307 

5 

0 

0 

0 

0 

0 

0 

D 

0 

0 

0 

6 

D 

0 

0 

0 

0 

3 

0 

0 

0 

4, 

7 

0 

D 

0 

0 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

0' 

1 

0 

0 

0 

o' ' 

0 

9 

D 

0 

a 

D 

D 

0 

0 

D 

13 

IT 

10 

D 

0 

0 

0 

D 

0 

0 

0 

0 

' 1 ' 

Total 

499 

38 

64 

l48 

1 — 

47 

131 

1 18 

1159 

38 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

iO - Undecided 
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Table 11. Mission 158 - Color film, filter 15 + 30B. Empirical Distribution 
Analysis with 2 equal frequency intervals. 


PREDICTED 


ACTUAL TYPE 


TYPE 

1 

2 

3 

4 

5 

6 

7 

8 

' 9 

Total 

' 1 

443 

T 

29 

50 

18 

72 

16 

70 

34 

739 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

56 

31 

35 

98 

29 

59 

2 

89 

4 

403 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 . 

6 

0 

0 

0 

0 

^ 0 

0 

, 0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

0 1 

0 

0 

0 

0 j 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

. 0 

10 

0 

0 

0 

0 

. 0 

0 

, 0 

0 

0 

0 

Total 

499 

38 

64 

148 

4T 

131 

18 

159 

38 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 12. Mission 158 - Color film, filter 15 + 30B. Empirical Distribution 
Analysis with 3 equal frequency intervals. 


ACTUAL TYPE 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 13. Mission 158 - Color film, filter 15 + 30B. Empirical Distribution 
Analysis with 4 equal frequency intervals. 


PREDICTED 

TYPE 


ACTUAL TYPE 


■9 


Total ' 




21 


55 


17 


66 


l6 


6k 


18 


707 


13 


l6 


11 


57 


Total 


1+99 


38 


61 + 


ll+8 


1+7 


131 


18 


159 


38 


2l* 


20 


l6 


28 


62 


21 


25 


6o 


2l+0 


13 


20 


53 


1+6 


10 


11 


10 


1 + 


LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 14. Mission 158 - Color film, filter 15 + 30B. Empirical Distribution 
Analysis with 5 equal frequency intervals. 


PREDICTED 

TYPE 


ACTUAL TYPE 


Total 


10 


452 


10 


22 


50 


15 


78 


l4 


64 


IT 


722 


0 


23 


17 


23 


65 


l4 


22 


31 


199 


0 


20 


Total 


499 


38 


64 


l48 


47 


131 


18 


159 


38 


16 


l6 


10 


i4 


33 


16 


54 


153 


15 


25 


LEGEND; 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 “ Water 

10 - Undecided 
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Table 15. Mission 158 -Color film, filter 15+ 30B. Empirical Distribution 
Analysis with 6 equal frequency intervals. 


PREDICTED 


ACTUAL TYPE 


X JIj 


C, 



. J 


( 

KJ 


XU ueu. 

1 

i+37 

7 

23 

42 

12 

73 

12 

55 

16 

677- 

2 

2 

2 

0 

0 



■ 0 

0 

0 

4 

3 

2 

0 

5 

0 

B 

2 

0 

4 

1 

l4- 

L 

- 1 

19 

12 

15 

46 


25 

B 

23 

5 

151- 

5 

7 

0 

1 

2 

Q 

5 

1 

B 

6 

0 

39 

6 

8 

1 

1 

5 

2 

10 

0 

6 

1 

34 

7 

0 

p 

1 

0 

0 

0 

0 

; 2 

0 

0 

2 

8 

19 

15 

18 

53 

13 

l6 

0 

6o 

-0 

194 

9 

0 

1 

0 

i 

, 0 

0 

0 

0 

2 

15 

18 - 

10 

5 

0 

1 1 

0 

0 

— 1 

0 1 

0 

3 

0 

9- 

Total 

k99 

38 

64 

l48. 



. 18 ,, 

159 . 

38 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 
k - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 16. Mission 158 - Color film,. filter 15 + 30B. Weighted Empirical 
Distribution Analysis with 4 equal frequency intervals. 

. _ f requencv/land use in cell 

Weig t - total f requency/land use 


PREDICTED 


ACTUAL TYPE 


TYPE 

1 

2 

3 

4 

5 

6 

7 

8 

9. 

Total 

1 

248 

■ 

3 

5 

2 

12 

3 

20 

2 

296 

2 

23 

l4 

B 

24 

2 

7 

0 

12 

0 

86 

3 

24 

9 

35 

49 

16 

22 

2 

53 

7 

217 

4 

16 

4 

0 

l6 

0 

3 

0 

7 

0 

46 

5- 

9 

5 

B 

13 

11 

l4 

0 

22 

0 

78 

6 

48 

3 

9 

21 

1 

36 

0 

17 

0 

135 

T 

86 

0 

4 

10 

13 

20 

11 

11 

5 

160 

8 

12 

0 

0 

2 

- 1 

2 

0 

9 

0 

26 

9 

30 

2 

5 

8 

1 

l4 

2 

8 

24 

94 ,, 

10 ■ 

3 

0 

0 

0 

0 

1 

0 

0 

0 

4 

1 

Total 

499 

38 

64 

l48 

47 

131 

18 

159 

38 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 
k - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 17. Mission 158 - Color film, filter 15 + 30B. Weighted Empirical 
Distribution Analys is with 4 equal frequency intervals. 

f requencv/land use in cell 
Weight - f requency/land use) ^ 


PREDICTED 


ACTUAL TYPE 



-L 

c. 



J 


I 


-7 

J-W 

1 

. 1 

396 

1 \ 

l6 j 

4o 1 

9 

Ll 

10 

k5 

12 

576 

1 

2 

17 

10 

2 

1 

16 

1 

6 

0 

11 

0 

63 

3 

10 

2 

"■ ' i 

l6 

2 

6 

11 

. 1 

1 

9 

5 

62 

U ^ 

13 

Ik 

18 

58 ■ 

Ik 

17 

1 

53 

0 

188 

5 

L 

0 

1 

7 

7 

9 

0 ! 

1 

13 

0 

Ul 

6 

26 

3 

8 

16 

1 

29 

0 

12 

0 

95 

T 

16 

0 

0 

k 

7 

' 9 

6 

3 

3 

48 

8 

12 

0 

0 

2 

1 

2 

' 0 

9 

0 

26 

9 

2 

2 

3 

3 

1 

6 

0 

1+ 

18 

39 

10 

3 

0 

1 

0 

0 

0 

1 

1 

1 

D 

0 

0 

1 . 

4 

Total 

k99 

i 

38 

6k 

, lU8 

1 

131 

, 18 

159 

38 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 
k - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 

10 - Undecided 
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Table 18, Mission 158 - Color film, filter 15 + 30B, Linear Discriminant 
Analys Is . 


PREDICTED 


ACTUAL TYPE 


TYPE 

1 

2 

L 3J 

4 ^ 

5 

6 

7 

8 

9 ' 

Total 

1 

118 

0 

■ 

■ 

4 

12 

0 

13 

1 

153 

2 

20 

3 

L 

5 

1 

6 

1 

B 

2 

46 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

k 

4o 

IT 

18 

54 

20 

35 

1 

44 

1 

230 

5 

1 

0 

0 

0 

1 

1 

0 

1 

0 

4 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

2k 

0 

1 

0 

1 

0 

1 

2 

0 

29, 

8 

k 

2 

6 

2 

1 

5 

0 

7 

B 

31. 

9 

59 

1 

5 

10 

5 

13 

9 

9 

21. 

132 

Total 

266 

23 

35 

75 

33 

72 

12 

80 

29 



LEGEND: 

1 - Forest 

2 - Crops 

3 - Plowed 

4 - Pasture 

5 - Idle 

6 - Abandoned 

7 - Orchard 

8 - Urban 

9 - Water 
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Table 19. Mission 158 - Color film, filter 15 +308. 
Discriminant Analysis. 


PREDICTED 

TYPE 


1*02 


97 


Total 


1+99 


ACTUAL TYPE 
2 

173 

1 + 70 ' 

61+3 


Total- 


575 


567 


LEGEND: 

1 - forest 

2 - nonforest 


Li near 
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Table 20. Mission 158 - Color film, filter 15 + 30B. 
Empirical Distribution Analysis with 3 equal frequency 
i ntervals . 


ACTUAL TYPE 


PREDICTED 




TYPE 

1 

2 . 

Total ■ 

1 

358 

136 

49 U 

2 

1^1 

.507 

648 

3 

, 0 

0 

0 





Total 

L99 ■ 

6U3 




LEGEND: 

1 - forest 

2 - nonforest 

3 - undecided 
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Table 21. Mission 158 - Color film, filter 15 + 30B. Linear Discriminant 
Analys is . 


PREDICTED 

TYPE 

1 

2 

1 

3 

4 

5 

1 

ACTUAL 

1 S ^ 

TYPE 
' 8 

9 

10 

11 

12 

13 

Total 

1 

46 

10 

21 

22 

0 

2 


, 0 

11 

0 

' 7 

2 

0 

123 

2 

11 

17 

10 

11 

3 

3 

2 

3 

B 

3 

8 

B 

0 

77 

3 

3 

B 

6 

B 

0 

- . 1 

1 

B 

5 

B 

2 

6 

0 

0 

42 

4 

0 

0 

0 

0 

i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

B 

B 

B 

5 

3 

1 

r ' ^ 

6 

2 

6 

0 

7 

B 

0 

4o. 

6 

0 

0 

0 

0 

0 

B 

° I 

0 

0 

1 

0 

0 

0 

2 

T 

7 

2 

8 

12 

14 

20 

50' 

19 

33 

1 

4i 

0 

0 

207 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11 


0 

1 

B 

2 

5 

3 

0 

B 

0 

B 

5 

1 

0 

28 

12 

2 

6 

7 

10 

1 

1 

■ 2 

8 

3 

B 

5 

5 

l6 

3' 

77 

13 

6 

5 

■ 

9 

0 

0 

0 

1 

B 

0 

2 

Bl 

0 ■ 

29 

Total 

80 

48 

i 

80 

23 

35 

75 

33 

. 72 

■ 

80 

26 

. 3 



LEGEND: 

1 - Pine 

2 - Pine hardwood 

3 - Bottomland hardwood 

4 - Upland hardwood 

5 - Crops 

6 - Plowed 

7 - Pasture 

8 - Idle 

9 - Abandoned 

10 - Orchard 

1 1 - Urban 

12 -Clear water 

13 •" Turbid water 
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Table 22. Mission 158 - Color film, filter 15 + 308. Empirical Distribution 
Analysis with 3 equal frequency intervals. 


ACTUAL TYPE 


PREDICTED 


TYPE 

1 

2 

! 

3 

4 

5 

1 1 

6 

7 

8 

9 

10 

11 

12 

L 

13 

] 

Total 

1 

11 4 

47 

42 

74 

■ 

6 

l6 

6 

32 

1 

5 

30 

■ 

2 

382 

2 

0 

B 

0 

0 

0 

0 

0 

0 

0 

, 0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

27 

23 

34 

58 

B 

10 

28 

6 

35 

7 

26 

6 

! 

1 

265 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

1 

0 

1 

0 

0 

0 

> 0 

0 

0 

0 

2 

7 

12 

3 

10 

l4 

23 

28 

75 

i6 

33 

1 

62 

3 

0 

280 

8 

0 

0 

. 0 

0 

0 

0 

0 

0 

0 

0 

0 

' 0 

0 

0 

9 

0 

8 

6 

4 

1 

2 

8 

7 

9 

2 

3 

2 

0 

52 

10 

0 

0 

0 

0 

0 

0 

1 

0 

' 0 

2 

0 

0 

0 

3 

11 

B 

0 

6 

9 

6 

17 

19 

12 

21 

1 

37 

9 

0 

l44 

12 

’ 0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

11 

0 

l4 

13 ■ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


; 0 

0 

1 ■ 

0 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

' 0 

B 

0 

Total 

i6o 

8l 

98 

i6o 

38 

64 

! 

l48 

i 

47 

131 

18 

159 

35 

i 3 



LEGEND: 


1 - P i ne 

2 - Pine hardwood 

3 - Bottomland hardwood 

4 - Upland hardwood , 

5 - Crops 

6 - Plowed 

7 - Pasture 

8 - Idle 

9 - Abandoned 

10 - Orchard 

1 1 - Urban 

12 - Clear water 

13 - Turbid water 

14 - Undecided 
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TREND AND SPREAD OF BARK BEETLE 


INFESTATIONS FROM 1968 THROUGH 1971 
R. C. Heller, K. A. Zealear, and T. H. Waite 

Bark beetles continue to concern forest managers throughout the 
United States; the most damaging insects are the southern pine beetle 
(Dendroctonus frontal is Zimm.) in the East and South, the mountain pine 
beetle (Dendroctonus ponderosa Hopk.) in the Rocky Mountain and Pacific 
Northwest States, and the western pine beetle (Dendroctonus brevicomjs 
LeConte) in the Southwestern States and California. The detection and 
location of infestations ("group kills" or spots) while they are small 
are still very important, so that control action can be started early 
before expensive epidemics get under way. For example, the costly 
mountain pine beetle outbreak on the Teton and Targhee National Forests 
and Grand Teton National Park during I 967 to 1970 (J_) might have been 
averted if over-mature trees had been removed in the early stages of the 
epidemic. 

Knowledge of the impact of the beetles on our national timber resource 
is still fragmentary and needs improvement. Studies such as those we have 
conducted in the Black Hills indicate that remote sensing improves our 
knowledge of the losses, where they are occurring, and does it more 
efficiently than methods used heretofore. 

By following the trend of how beetle populations operate over a 
time period, we can expect to learn where new infestations are likely to 
occur and to predict roughly the size and location of new infestations. 
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For example, small-scale aerial color photographs (1:32,000) permit the 
forest manager to evaluate current losses more accurately over larger 
areas than he could do by ground reconnaissance methods. When the small- 
scale photos are used with medium-scale photography (1:8,000) and ground 
examination, we have a most efficient survey design { 2 ^. One of the 
by-products of the NASA-Forest Service cooperative studies is the collec- 
tion of four years of bark beetle infestation data over one relatively . 
small area (1.6 by 5 kilometers) where little or no effort was made to 
control the beetle epidemic. 

The study area in the Black Hills was selected because the timber 
type, is relatively uniform -- more than 80% ponderosa pine. The remain- 
ing species are white spruce (Picea qlauca (Moench) Voss.), trembling 
aspen (Populus tremuloides Michx.), and white birch ( Betula. papyri fera 
Marsh). Also, there is only one generation of beetles per year in this 
area, which causes only one new generation of dying and discoloring 
trees to be present at any one time. This "once-a-year" change simplifies 
the separation of "old group kills" from "new group kills" on aerial 
photographs. From August through September of each year, these newly 
discolored trees offer relatively large yellow-red targets which have 
a fairly high contrast ratio to the uniformly green healthy forest. 

Since the infestations occur in random sizes and locations, they become 
ideal resolution targets which we have used to answer questions on 
film-filter-scale combinations. These known targets are also expected 
to be used for possible detection by the Earth Resources Technology 
Satel 1 i te. 
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PROCEDURES 


Aerial color photography was taken for three purposes over selected 
areas in the northern Black Hills on August 11 and 12, 1971. First, we 
wanted to rephotograph the 1.6 by 5 kilometer study area mentioned 
earlier for following the trend of the beetle epidemic. Secondly, we 
wished to locate prospective ERTS sites in the event we were accepted to 
undertake our proposed ERTS research. Finally, we wished to establish 
the usefulness of aerial photography to identify levels of tree density 
which would not support beetle infestations. These two latter studies 
will be reported elsewhere in this progress report. 

Only one scale of photography was flown -- 1:8,000 — with a Zeiss RMK 
21/23 from the Forest Service Aero Commander. Both medium-speed Ektachrome 
(Type 2448) and Aerocolor (Type 2445) were exposed over most sites. The 
medium-speed Ektachrome was processed in our own Berkeley facility to 
a positive transparency. The color negative film was sent out for develop- 
ment in a continuous processor. Color prints were made on a color Loge- 
tronic printer at the Forest Service color lab in San Francisco. 

All i nterpretation was carried out on the color transparencies; 
resolution and color fidelity were better on the transparencies than on 
the prints. All newly discolored trees, which were yellow to yellow-red 
i n Augus t, were circled in ink directly on the transparencies. The 
groups were numbered and the trees counted in each group. Old-killed 
trees (1970 faders) had to be distinguished from the new faders (discolored 
in 1970 srid this, too, could be done more accurately on the transparencies 
than on the prints (Figure 13). Finally, within the 1.6 by 5 kilometer study 
area, all data on the transparencies were transferred to the color prints 
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Figure 13. A portion of the 1.6 by 5 km. study area near 
Lead, S. D., is shown in this color print made from Aerocolor 
film. Scale 1:8,000, Zeiss RMK 21/23 camera. The newly dying 
trees are a yellow hue and are circled on this print. The 
infestation number shown near the circled trees can be checked 
in Table 8 for ground tree count. Note hue discrimination of 
older killed trees (yellow-red to brown) from healthy (green 
to green-yellow). 
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BLACK HILLS BEETLE STUDY 
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Ground Tally 


Date 9/9 / 7 / 


Crew 





Sheet of 

Figure 1 A. Sample of field form used to measure size of each infestation 
and tally of trees within each infestation. 
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for field use. The Interpreters were instructed to include all possible 
discolored trees so that there would be little likelihood of omission errors 
occurring. A total of 390 infestations were located on the photos. 

The color prints were used stereoscopical ly in the field and were 
of incalculalable benefit in navigating through the timber. All 390 
suspected infestations were ground checked for (1) presence of beetle 
attack, (2) number of trees within each group, (3) bole diameter in inches, 

Q 

(4) crown class, and (5) the size of the infestation in meters (Figure l4) . 

RESULTS 

All trees which were found to be beetle infested were plotted on an 
overlay map of the study area (Figure 15). Table23 lists these infestations 
by numbers of trees per group and by longest dimension of the infestation 
in meters. While 390 infestations were located on the photos, only 253 
were found to be 1971“faded infestations when checked on the ground. 

Most of the 137 commission errors (390 - 253 = 137) were found to be 
single trees which probably faded late in 1970. The western section of 
this map, where the greatest amount of tree killing has occurred during 
the past four years (1968 through 1971). is shown in overlay form (Figure 
16) for each of the years. Note particularly how the larger infestations 
have enlarged from one small location to merge with other large infesta- 
tions. Hundreds of large trees, 25 to 30 meters tall, which were killed 
in 1968 are now down on the ground, changing the forest composition and 
environment for these sites. 

Table 24 is a listing similar to Table23 but includes only the infesta- 
tions shown on the top overlay of Figure 16. The small numbers shown next 

^Crown class is a silvicultural term relating the position of a tree 
to its neighbors. For example, the crown of a dominant tree receives 
direct sunlight from above and from three or more sides and is easily seen 
on air photographs. The crown of an intermediate tree receives only light 
from above and is often missed during photo interpretation. 
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Table 23 - List of 1971 mountain pine beetle infestations on study area 2 by number of trees per infestation and 
longest dimension of the infested area. 
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Figure 15. Aerial mosaic on left of 1.6 by 5 km. study area near 
Ledd, S. D. On right, 253 infestations of various sizes (3 to 20 
meters) are plotted at same scale. Outlined area on right refers 
to photo coverage shown in Figure 13 and the overlay shown in 
Figure 16. 




Figure 16. Overlay of four years of mountain pine beetle infestations 
within selected portion of study area (same coverage as photo shown 
in Figure 13 ). Hardwoods are shown on base map as being crossed- 
hatched. 1971"killed trees are shown on top overlay in red; 1970 in 
orange; 1969 in blue; 1968 in green. Note the grouping of the infes- 
tations. Most of the pine timber has been killed in this small area. 
Size and number of each infestation in 1971 are referenced in Table 2k. 
Scale 1:8,000. 
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Table 24. Infestations by spot number, number of trees, and longest 
dimension of the infestation. 


Spot 

No. 

No. of 
Trees 

Length 

Ft. 

m. 

Spot 

No. 

No. of 
Trees 

Length 

Ft. 

m. 

252 

29 

150 

45 

354 

236 

470 

143 

316 

948 

550 

168 

355 

1 

10 

3 

317 

7 

30 

9 

357 

24 

132 

41 

323 

9 ' 

60 

18 

359 

81 

210 

64 

324 

10 

90 

27 

360 

1 

10 

3 

325 

3 

30 

9 

361 

14 

40 

13 

cr\ 

CM 

cc\ 

2 

15 

5 

362 

7 

32 

10 

330 

1 

10 

3 

363 

337 

470 

173 

331 

61 

200 

61 

364 

1 

10 

3 

334 

90 

265 

81 

365 

347 

950 

290 

335 

1 

10 

3 

368 

33 


27 

336 

3 

30 

9 

369 

104 

330 

9 

344 

7 

55 

17 

371 

10 

70 

21 

345 

1 

10 

3 

372 

39 

200 


346 

1 

10 

3 

373 

45 

205 

62 

347 

1 

10 

3 

374 

1 

10 

3 

00 

1 

10 

3 

376 

303 

375 

114 

349 

4 

30 

9 

388 

19 

190 

58 

350 

1 

10 

3 

ON 

00 

2 

25 

8 

351 

1 

10 

3 

390 

2 

15 

5 

352 

8 

80 

24 
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NUMBER OF TREES PER INFESTATION - P.l. COUNT 



Figure 17- Regression of ground counts to photo counts. Confidence bands 
at the 95% level are shown on each s i de of regression line. For 1;8,000 
scale color transparencies. Includes all crown classes. 
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to the circled infestations on the color print (Figure 13) refer to 
the listed infestation numbers in Table 24, 

Photo interpreters cannot see all trees in a stand of timber. 

Usually fewer trees are counted on photos as scale becomes smaller. 

Also, sma 1 1 -crowned single trees are frequently missed both on the 
photos and on the ground. These facts have been documented in several 
other studies (3,, 4), However, usually a relationship between ground 
counts and photo counts can be developed from which a ground count can 
be estimated and a level of confidence can be expressed for any selected 
photo count. A linear regression was developed for these data (Figure 17) 
which should be helpful in making total mortality estimates over large 
areas when 1:8,000 color photography is used as the sampling frame. One 
can use Figurel7by entering the ordinate (or photo count) reading across 
to the regression line and down to the ground count. For example, a photo 
count of 25 trees per infestation would yield a ground count of 38 trees 
with an estimated possible error of + 6 trees, 19 out of 20 times. The 
reason for the increase of the ground count over the photo count is that 
the photo interpreter cannot see all trees in the overtopped and inter- 
mediate crown classes. 

Figure l8 is a plot of three regressions. The lowest curve (all 
crown classes) is identical to Figure 17 on which the confidence bands 
have been superimposed. Note that as the suppressed and intermediate 
tree crown classes are removed from the ground count, the regression 
slope approaches 1.0. Undoubtedly the confidence limits would be even 
closer to the plotted regression than in Figure l7. 
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Figure 18. Comparison of ground tree counts to photo tree 
counts of discolored ponderosa pine infestations on 1:8,000 
color transparencies. The lower line is identical with Figure 
17 and includes all crown classes. The middle line -- dashed 
includes dominant, codominant, and intermediate crown classes. 
The upper solid line is for the largest trees in the dominant 
and codominant crown classes only. 
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When assessing the detectability of stressed trees from ERTS imagery, 
the map of 1971 infestations will be extremely helpful in providing 
natural resolution targets whose size we know quite accurately. 

These data will be useful until July 1972; however, by mid-August 1972, 
new infestations will appear and provide a different set of resolution 
targets. These new targets will have to be assessed on the ground at 
that time. Additional discussion on the probability of detection of stressed 
trees from ERTS imagery can be found in the report by Weber and Heller, 
wh i ch fol lows . 
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AERIAL PHOTOS DETERMINE OPTIMUM LEVELS 


OF STAND DENSITY TO REDUCE BARK BEETLE INFESTATION 

by 

T. H. Waite, R. C. Heller, and R. E. Stevens 
INTRODUCTION 

Large bark beetle infestations have been and will continue to be 
costly until adequate management practices are adopted. They are costly 
not only in terms of timber production but are damaging to watershed, 
wildlife, and recreational resources as well. 

For years experienced field men have noted an apparent relationship 
between low stand basal area and resistance to bark beetle attack. Field 
observation of ponderosa pine ( Pinus ponderosa Laws.) stands in the Black 
Hills indicates that stands which have been thinned appear resistant to 
beetle attack. This same relationship was evident on our aerial photo- 
graphs of the area. 

Because aerial photographic techniques may permit more rapid selection 
of stands needing thinning, this study is designed to show the relationship 
between beetle attack and stand basal area. It should provide evidence to 
determine if thinning is a viable treatment for bark beetle control. 

The objectives of this study were twofold. The first was to determine 
whether aerial photographs can be used to identify various stand density 
levels. The second was to determine what density levels will reduce the 
likelihood of beetle attack. 
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PROCEDURES 


In August 1971. new photography was flown over the Black Hills study 
area (1:8,000 Aerocolor 2445 and MS Ektachrome 2448). Five separate areas 
were flown in order to insure broad coverage of infested, uninfested, and 
thinned areas. 

To obtain a representative sample of the study areas, two photographs 
were selected at random from each of the five areas photographed. In 
study area II an exception was made. The photos were selected purposely 
to include areas of heavy infestation. A grid having 63 one-inch squares 
was placed over the selected photos. The outside 1- x 9- i nch strip on 
each side of the photograph, and parallel to the flight line, was elimi- 
nated in order to reduce distortion in the x direction. Each square inch 
(approximately 10 acres) of the sample was stratified into one of two 
categories: (1) pine type with no infestations and (2) pine type with 

infestations; nonpine type was not sampled. Each stratum was listed by 
the photo of origin and the identifying square inch. From this list, 

40 sample one-inch squares were chosen at random from each stratum. These 
one-inch squares were further divided into four 1/4-inch (2 1/2-acre) squares 
one was selected at random and listed by corner (NE, SE, SW, NW) , square, 
and photo number. The final random sample of twenty 1/4-inch squares 
was drawn from each stratum, giving a total field sample of 40 (two were 
rejected in the field as being out of type, giving a field sample of 38; 

22 were infested, and 16 were uninfested). 

The sample of the thinned stands was obtained in a different manner. 
Again, all 1971 photography taken of each of the five areas was considered. 

The thinned areas on the photos were delineated. Recently thinned areas 
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(defined by the presence of fresh haul roads, equipment or active landings) 
were not considered. On each of the thinned areas two 1-acre plots were 
selected at random (avoiding road edges and openings). A pinprick was 
placed on the photo In the center of the plot, and estimates of crown 
closure (density) and average crown diameter were made using a crown 
closure comparator (J_) . In all, sixteen plots were located on thinned 
areas. 

All 5^ plots (38 nonthinned and 16 thinned stands) were visited on 
the ground. The stereo color prints used in the field were of incalculable 
value in locating the points. Each of the 54 plots was recorded in the 
same manner. The center point pricked on the photo was visited, and the 
photo and point numbers were recorded on the field form. A ten-point 
plot cruise using a 40 B.A, prism was done at each sample point. In 
addition, at each point the the crown diameter, crown class, and 

species were recorded; species other than pine were included ip the tree 
count only, as it was felt that these factors contributed to basal area 
and thus to stress. At 5 of the 10 points, dominant tree height was 
recorded . 

RESULTS 

Three separate analyses were made of the data. 

1. To determine the relationship of crown comparator measure- 
ments (crown closure and crown diameter) to basal area. Are they highly 
correlated? To what accuracy? 

2. To test the null hypothesis that thinned stands are from the 
same population as natural stands. 
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3. To test the null hypothesis that the basal area of non- 
infested stands is from the same population as the basal area of in- 
fested stands. 

In order to determine the relationship between crown closure and 
crown diameter measurements on the photographs and stand basal area 
measurements on the ground, two simple linear regressions were made 
(F i gures 1 9 and 2D) . F i gure 19 shows the linear relationship between basal 
area and crown closure. The regression coefficients were tested using 
Student's "t" test and were found to be significant at the 95% confidence 
level. This means that if there is in fact no relationship between the 
two variables, given this many observations, the chance of getting co- 
efficients this different from zero is 5%. A more intuitive explanation 
would be: This result indicates that a si gnificant amount of the vari- 

ation in the basal area values is explained by the fitted linear regres- 
sion. However, its value as a predictive tool at this point is questionable. 
For instance, the coefficient of determination r^ was only 0.47, Although 
this value is significant at the 95% confidence level, it implies that 
only 47% of the variation in the basal area values is associated with 
the crown comparator measurements. Furthermore, the standard error of 
the basal area estimate is +40.72 square feet. This value is quite high 
and is an indication of the rather large variation of the basal area 
values about the mean. 

F i gure 20 shows the linear relationship between basal area and crown 
diameter. As we expected, the regression slope coefficient is negative, 
indicating that as stand basal area decreases crown diameters increase. 

Both regression coefficients were significant at the 95% confidence level. 
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PLOT BASAL AREA 



Figure 19. A simple linear regression showing the relationships of stand 
basal area as measured on the ground to average crown cover as measured by 
the crown comparator. 
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Figure 20. A simple linear regression showing relationship of stand basal 
area as measured on the ground to average crown diameter as measured by the 
crown comparator. 
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The Implications in terms of the model are the same as that for Figure 19> 
i.e,, a significant amount of variation in basal area can be explained 
by the fitted linear regression. The coefficient of determination, r^, 
for these data is only 0.29 — indicating that only 29% of the variation 
in the basal area values is associated with the crown diameter measure- 
ments. However, this value is significant at the 95% confidence level. 

On this basis we reject the null hypothesis that there is no relation- 
ship between the variables. In addition to the low r value, the large 
variation about the fitted line is further substantiated by the high 
standard error of the basal area estimate which is +49.27 square feet. 

Figure 21 is a plot showing the relationship between basal area and 
d.b.h. for both infested and uninfested plots. All thinned plots were 
uninfested and as such were included in the uninfested group. ^ In view 
of the distribution of the plotted data, it was decided that a .simple 
linear regression analysis was inappropriate. A Student^s “t” test was 
done on the differences between mean d.b.h. *s and mean basal areas of the 
infested and uninfested plots. The mean d.b.h. for the infested plots is 
12.64 inches and for the uninfested plots is 9.63 inches. Student's "t" 
test indicates that this difference is significant at the 95% level. 

This means that if there is no difference between . these groups, given 
this number of observations, a difference this large would arise 5% of 
the time. In other words, the probability is very high that these groups 
are from different populations. The mean basal area for the infested plots 
is 149.33 square feet and for the uninfested plots is 113.64 square feet. 
Again, Student's "t" test indicates that the difference is significant 
at the 95% confidence level. The implication is the same as that for 
the mean d.b.h. 's, i.e., the probability is very high that these two 
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Figure 21. A plot showing the distribution of infested, uninfested and 
thinned stands by average basal area and average stand diameter. 
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groups are from different populations. 

CONCLUSIONS AND DISCUSSION 

The crown comparator has shown itself to be a photo interpretation 
topi with substantial development potential. Further experience in its 
use should reduce the variability found in this study. The elementary 
analysis conducted on the thinning data indicates that reducing stand 
basal area will increase resistance to beetle attack. However, this 
phenomenon is probably more complex than our elementary analysis indicates. 
Our studies in the future will look toward a more complex analysis of the 
effects of basal area on insect attack. For instance, a multivariate 
approach would be suitable for this purpose. It would allow us to 
determine the contributions made by other variables such as slope, aspect. 
Soil, and micrpsite to the susceptibility of stands of beetle attack. 

Some of these data are available to us, and we expect to continue the 
a na 1 ys i s . 
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THE USE OF AIRBORNE SPECTROMETERS 


AND MULTISPECTRAL SCANNERS FOR PREVISUAL DETECTION 
OF PONDEROSA PINE TREES UNDER STRESS FROM INSECTS AND DISEASES 

F. P. Weber 

Several times during the last seven years we have looked at the 
application of opt ica 1 -mechani ca 1 infrared line scanners for previsual 
detection of stress and loss of vigor in pine trees. In the original 
feasibility study undertaken by the Forest Service ( 4 ), we tried fifteen 
individual spectral channels of data from the University of Michigan's 
original AAS -5 scanners. Because the focus of the supporting field study 
was to detect temperature differences between healthy and low-vigor trees, 
a group of forty-year-old red pine trees ( Pi nus res i nosa Ait.) were arti- 
ficially stressed with injections of an herbicide. A temperature rise in 
the affected trees was measured in the field and was related to the 
decline in tree vigor. The same temperature difference was recorded 
by the Michigan thermal scanner, prior to the time the stressed trees 
were detected by a change in color on large-scale color and color infra- 
red photography. 

Because of the success with the feasibility study in Michigan, we 
established a cooperative pilot study with NASA in the Black Hills in 
1966 to determine if trees attacked by the mountain pine beetle (Dendroc- 
•tonus ponderosae Hopk.) and an associated blue stain fungus could be 
detected on thermal infrared imagery prior to identification on large- 
scale color photography. Test sites were established and trees were 
instrumented to measure data on the biophysical responses and physio- 
logical differences of healthy and beetle-attacked ponderosa pine (^) , 
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The airborne f i re . detect i on and mapping group from the Northern 
Forest Fire Laboratory in Missoula, Montana, flew our Black Hills test 
sites in 1966 and 1967 with their bispectral thermal infrared scanner. 
Missions were flown both in daytime and at night and with various meteoro- 
logical conditions. In spite of the concerted effort by our joint research 
groups, we were not able to identify the beetle-attacked trees on thermal 
infrared imagery (J_, 2 ), 

In 1968 and 1969, we again worked closely with the Infrared and 
Optics Laboratory (IROL) of the University of Michigan. They had developed 
an improved scanner capability with the addition of a 12-channeI spectro- 
meter. This development opened the door to sophisticated mu 1 1 i spectra 1 
processing, at least for data in the visible and near infrared portions of 
the spectrum. New test sites were established in the Black Hills to provide 
base-line biophysical data and calibration to support the Michigan multi- 
spectral scanner data. 

Although our results with Michigan demonstrated the capability to 
identify bark beetle infestations, two findings were most noteworthy. 

First, although thermal infrared (8 to 14 ym) data contributes to the 
detection of beetle-attacked pine, reliable detection requires the use 
of s i mu la taneous visible, reflective infrared, and thermal infrared data 
collected in narrow-band channels selected from a broad-band mul tispectral 
system (0.4 to 14.0 ym). Secondly, hybrid mu 1 1 i s pectra 1 processing is 
required which utilizes the best features of a di gi tal computer for 
optimum channel selection in concert with a specialized analog computer 
which actually does the recognition process (J) • 
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Two recent instrument developments improve the potential for furthering 
our stress detection capabilities in 1972. Both the Bendix 24-channel multi- 
spectral scanner and the Michigan (M-7) 12-channel mul tispectral scanner 
incorporate design features which we have determined are necessary to 
Improve our present level of success in previsual stress detection. The 
principal improvement in both systems is the integrated design which 
permits simultaneous registration of all spectral channels of data. This 
we know is the single most important requirement of our work. 

Until now our position has been that for optimum detection of forest 
stress conditions where automated feature recognition is desired, it may 
be necessary to design and construct special systems. There have been few 
alternatives to this belief until now. The requirement of high spatial 
resolution for coniferous trees with small crown diameters as well as for 
deciduous trees with large crown diameters has been achieved with high 
sampling rates. Moreover, this workwoujd have to be done wi thout sacr i - 
ficing s ignal-to-noise ratio (S/N) capability which relates to the minimum 
detectable differences in temperature. In general, improving the spatial 
resolution for a scanner by a factor of 2 reduces the detector sensitivity 
by a factor of 4. This reduction means that a 3 milliradian (mr) resolu- 
tion system (with, for example, a 0.5^F Noise Equivalent Temperature sensi- 
tivity (NET), when converted to 1.5 mr systems by changing detectors may 
have only a 2^F NET. As has been shown by extensive ground measurements, 
this temperature is within the normal range of differences between .stressed 
and unstressed trees. Therefore, such a conversion would not provide 
sufficient S/N capability for reliable identification of tree conditions. 

This resolution dilemma can be alleviated now i n two ways. First, 
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we can use one or more of the new systems which provide truly registered 
mul t ispectra 1 data in narrow spectral bands over a broad portion of the 
spectrum. Secondly, we can specify flights at low altitudes to achieve 
better ground resolution. 

As part of our continued participation in the Earth Resources Program 
during 1972-73, we anticipate getting several supporting flights by the 
NASA C-I 30 aircraft with the Bendix 24-channel Multispectral Scanner Data 
System (MSDS) . Our program is in part designed to test the capability of 
the MSDS for detecting green beet le- i nfes ted trees in the Black Hills. As 
a minimum, we expect to demonstrate large gains in the area of automatic 
thematic mapping and landscape feature recognition. Additionally, the 
MSDS flights will be a coordinated part of the ERTS experiment for the 
Black Hills test site (226A) . Baseline biophysical and calibration data 
needed for processing the MSDS imagery will be collected at the test site 
and transmitted to ERTS for correlation with satellite-collected imagery. 

The focal point of the MSDS missions in early 1972 (within the 
3,000-square-mile ERTS test site) is the site selected for establishment 
of the Data Collection Platforms (DCP) . This location (n44^ 16‘, 103^ 

47*W) was selected primarily because it includes new bark beetle activity 
of the type required for a good test of the 24-channel scanner. The 
general location of the new study area is shown in Figure 22. By cross 
checking with Figure 26(page 111 ) of the section of this report by Weber and 
Heller, the reader can see how the MSDS study site and DCP location tie 
in with a larger portion of the ERTS experiment test site.. 

Within the MSDS study site, we have located one group of 54 large 
ponderosa pine trees which were attacked by the mountain pine beetle in 
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Figure 22. The location of the MSDS study site and the data collection 
platforms (within the larger ERTS test site) are shown on this 1:55,000 
scale photography. 
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late August 1971. This is in an area where there have been few beetle 
infestations for several years; this makes it ideal for the test. The 
newly attacked trees were green when the site was selected in September 
1971 and will remain that way until after the first requested flight 
support (April 1972) by the NASA C-130 aircraft. Subsequent support 
flights throughout early spring and summer should provide a good test 
of the capability of the MSDS for previsual detection of stress in 
conifers. 

An important part of the pre-experiment, study site establishment 
was the construction of a three-tower tramway system in September 1971 
(Figure 23). The tramway was established so that the south half of the 
north-south sample line runs, above healthy pine and the north half above 
the newly attacked pine. The purpose of the tramway system is to provide 
a means of transporting sensors for measuri ng spectra 1 radiance and other 
important biophysical data. Sensor data will be recorded much the same 
as that reported by Weber and Wear (b) . In addition, some data that 
are recorded by our digital data logger will be multiplexed and trans- 
mitted via two OCR's to ERTS . 

In past studies we have measured many radiation parameters to discover 
the implications to remote sensing of reflected and emitted energy differences 
between healthy and stressed trees. Always, a missing link has been the 
lack of spectral energy differences measured at the target. This past 
season we have built and calibrated a small portable field spectrometer 
which currently operates with four ERTS-matched channels. The first 
prototype model was constructed with the intent that it be suspended 
from the tramway system to provide baseline spectral information both 
for ERTS and for the MSDS flights. The operating prototype of the 
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Figure 23. The alignment of the three-tower tramway system (A) and the 
DCP location (B) are shown in this September 1971 aerial photo. The group 
of 5^ new beetle-attacked trees is also shown (C) . 
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field spectrometer was checked out in the laboratory and the results 
were better than expected. It will be field tested the first week in 
November as a coordinated effort with the flight of the Michigan M-7 
scanner system over our Atlanta test site (NASA Site 217)- 

The 1:8,000 scale aerial color photograph (Figure 24), taken of the 
Black Hills MSDS study site in August 1971> shows several of the landscape 
features which will be instrumented for the coordinated ERTS-MSDS experi- 
ment. Spectral data, in addition to other critical biophysical and 
physiological information, will be gathered at a range site and a 
soil outcrop site in addition to the forest site. The range site was 
selected by R. S. Driscoll, Principal Range Ecologist from the Rocky 
Mountain Forest and Range Experimental Station, who will assist In collec- 
tion and analysis of range data, 

We look forward to the opportunity In 1972 for a field trial of the 
Multlspectral Scanner Data System for previsual detection of green bark 
beetle-attacked trees In the Black Hills. We have spent much of the last 
seven years defining our multlspectral scanner requirements and are eager 
to try the Bendix 24-channel system and the Michigan 12-channel M-7 system. 
During the same period we have developed a good scheme for evaluating bark 
beetle Impact on forests, using a combination of high-resolutlon aerial 
color photography and probability statistics (2). However, current 
photographic analysis techniques do not lend themselves readily to auto- 
matic data processing, Although some automatic photo analysis can be done 
with film, It can be done more easily If the data are stored Initially on 
magnetic tape. Moreover, multlspectral scanner techniques seem to provide 
the best means for previsual detection of stress In forests. The efficiency 
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Figure 2k. The focal point of the Black Hills MSDS study site is shown 
on a 1:8,000 scale aerial color photograph taken in August 1971. The 
pasture in the center of the photo and the light, soil-rock outcrop just 
above the field will be instrumented in addition to the forest site Just 
to the right of the upper end of the pasture. 
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and reliability of scanner techniques in comparison to photography 
remains to be demonstrated. 
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ESTABLISHMENT OF EARTH RESOURCES TECHNOLOGY SATELLITE 


TEST SITE IN THE BLACK HILLS AND IDENTIFICATION OF NATURAL RESOLUTION TARGETS 

F. P, Weber and R. C. Heller 

The Black Hills National Forest test site (NASA site 149) encompasses 
a 3,000-square-mile forest and rangeland area in western South Dakota and 
eastern Wyoming (Figure25). For the past seven years we have maintained an 
intensive aerial and ground study in the Black Hills for learning how to 
use mul tispectral sensing techniques to detect trees under stress from 
insect or disease attack. During that period of time we have observed 
several phenomena of change, such as the manifestation of dead trees from 
insect attack, that have occurred over a large enough area or resulted 
in sufficient spectral change to suggest the feasibility of change detec- 
tion from space. 

In mid-September 1971, after an extensive survey of the northern 
Black Hills, we selected a location for the ERTS data collection platforms 
(DCP) . The site -- n 44° 16', 103° 47'W -- is a central location which 
embraces all of the physical characteristics of the landscape we hope 
to study with ERTS imagery. Features such as range grass, soil-rock out- 
croppings, and beetle-attacked and health pines are close enough to the 
DCP site so that biophysical data from each target, e.g., spectral 
radiance, can be measured and recorded via land lines to the central 
data collection site. 

An important consideration in the site selection was the location of 
a new mountain pine beetle infestation in the area to be used as an indi- 
cator of tree stress. The new infestation includes 54 large trees grouped 
in roughly a circular pattern with an area of approximately 1,260 square 
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BLACK HILLS, SOUTH DAKOTA 
TEST SITE 149 
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Figure 25. The Black Hills National Forest test site in western South 
Dakota and eastern Wyoming. Notice the location of the ERTS data col- 
lection platform (DCP). 


106 


camera lens (j_) . Each distribution was integrated by the trapezoidal 
rule from 400 to 700 nm and normalized to the same bandpass value. 

The def i ni tion is 


( 1 ) 


where 


680 

j [D400 + O700I + 2 0. = 25.00 

X=420 ^ 

D, = log relative daylight irradiance, 
A 10 


X = wavelength in nanometers. 


Another way of expressing the meaning of these data is that they represent 
the relative daylight distribution passing through a lens onto film, normalized 
to the average value in the 400-700 nm range. 

For infrared-sensitive film which is almost always used with a minus- 
blue filter, we attenuated the curves by the density of the Eastman Wratten 
12 filter and chose 500-920 nm as the representative bandwidth. The normal- 
ization equation is 

900 

(2) j [D 500 + D 920 ] + Z ^ = 50.00. 

^=520 

The data of reference 5 extend only to 83 O nm so that distribution was 
extrapolated to 920 nm using the sunlight data of reference 2. Skylight 
does not appreciably alter the shape of the daylight curve between 830 and 
920 nm. 

Each daylight distribution D is denoted by the reference number from 
which the data were drawn. DIA is a composite sun + sky distribution as 
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meters. The beetle attacks occurred in late August 1971, so that at 
the time of the expected launch and first data collection by ERTS, the 
attacked trees will still be green. However, during the first five 
months of the ERTS experiment, the study trees should discolor from 
green to yellow-green and finally to yellow-red by August 1972. 

During our September visit to the Black Hills, 10 large groups of 
beetle-killed trees were selected within our existing study area II for 
observation and study during the ERTS experiment (Table 25). See Figure 
13in the paper by Heller, Zealear, and Waite of this report showing these 
infestations (page 7^). After we had carefully examined the infestation spots 
on aerial color photography, each spot was assigned a relative probability 
of detection on ERTS imagery. The assignment of a relative value -- low, 
moderate, or high — is based on three characteristics: (1) spot size, 

(2) spot geometry, and (3) spectral contrast of the spot in relation to 
the surrounding landscape features. 

The assignment of detection probabilities based on appearance on the 
August 1971 photography is a relative entity attached to a dynamic process. 
That is, the trees continue to change in color and texture between August 
1971 snd launch time, March 1972. The prediction considers that beetle- 
killed trees that were yellow-red this past summer may be indistinguishable 
from other old-killed trees by the spring and early summer of 1972. 

Spots in this situation would have a low detection probability in contrast 
to spots which were completely surrounded by healthy trees. Our main 
purpose here is to determine various beetle spot sizes that can be recog- 
nized on space imagery relative to spectral contrast with the surrounding 
landscape features. This objective implies that one moderate-size spot in 
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Table 25. Numbers of trees and size of mountain pine beetle infestation 
spots identified for the ERTS experiment. 


Spot-'^ 

Number 

Number of 
Attacked 
Trees 

Longest 

Dimension 

(meters) 

Area 

(meters^) 

Probability^' 
of Detection 

176 

108 

38 

1,090 

low 

229 

214 

214 

9,750 

moderate 

309 

173 

124 

4,840 

low 

311 

143 

76 

2,240 

low 

316 

948 

225 

10,125 

moderate 

354 

236 

147 

11,700 

high 

363 

337 

147 

5,900 

moderate 

365 

347 

290 

17,700 

high 

369 

104 

no 

3,380 

low 

376 

303 

155 

6,975 

moderate 


]_/ These infestations ("group kills") are circled and numbered on the color 
pri nt (Fi gure 1 3) . 

"U Probability of detection is based on total size and shape of beetle spot 
in addition to spectral contrast with healthy pine, old-killed pine, and'other 
landscape features. 
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a green forest of healthy pines may be more readily identifiable than, 
say, a group of large spots in an area of intense insect activity where 
spot aggregation is prevalent. 

Initially, there are three additional physical features which we 
would like to be able to identify and evaluate on ERTS imagery in addition 
to bark beetle-killed trees. They are: (1) range grass and pasture land 

(both natural and planted), (2) soil and rock outcroppings, and (3) water 
storage areas and water courses (Figure 26). Within our ERTS study site 
we have identified three natural targets within each of the three cate- 
gories (Table 26). Once again, each target was assigned a detection proba- 
bility based on size, shape, and spectral contrast with surrounding features. 
The evaluation for deriving the detection probabilities was done on 1:220,000 
scale color photography obtained from the August 1969 RB-57 flight over the 
Black Hills. 

Scientists worldwide look forward to the availability of the first 
satellite-collected space imagery next year. We have tried to establish 
a forest test site, based on our experience in the Earth Resources Program 
in the Black Hills, which will reveal any possible utility to resource 
managers from current space imagery. We feel competent to make the evalu- 
ation of applicability; however, it remains to be seen just how good and 
useful the ERTS imagery will be. 


no 


Figure 26. A portion of the Black Hills test site for the ERTS experiment 
is shown on 1:110,000 RB-57 photography. Image annotations show several of 
the targets which are listed in Table 26. These targets represent the three 
major landscape features of interest, other than bark beetle-killed trees, 
which will be our focus for ERTS image analysis. 
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Table 26 

• Black Hills landscape 

features selected as 

natoral resolution 

targets 

for the ERTS experiment. 




Longest 


Probability^^ 

Natural 

Dimension 

Area 

Target 

(meters) 

(meters^) 

of Detection 


Range grass 
Pasture A 


(planted) 

403 

64,722 

moderate 

Pasture B 
(planted) 

905 

200,093 

high 

Pasture C 
(natural) 

2,262 

135,720 

low 

Soil-rock Outcrop 




Gravel pit 

295 

68,349 

moderate 

Mine tailings 

940 

197,400 

high 

Rock outcrop 

195 

11,700 

low 

Water 




Reservoir 

8,850 

1,371,600 

high 

Lake 

440 

152,053 

moderate 

River 

7,280 

- - - - 

low 

J_/ Probabi 1 i ty 
i n addi tion to 

of detection Is based 
spectral contrast with 

on tota 1 si: 
surround! H9 

ze and shape of target 
features . 
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DETECTION OF ROOT DISEASE IMPACTS ON FOREST STANDS 


BY SEQ.UENTIAL ORBITAL AND SUBORBITAL MULTIS PECTRAL PHOTOGRAPHY 

John F. Wear 
INTRODUCTION 

The feasibility of detecting and measuring the impact of root-rot 
diseases on our forest resources by applying remote sensing techniques 
from aircraft and earth satellites is an important study of the Remote 
Sensing Work Unit in the U. S, Forest Service. Earth resources managers 
are greatly concerned about the continuing timber losses over extensive 
forest areas caused by forest diseases. Satisfactory aerial surveying 
techniques are needed to acquire data on the distribution of disease 
centers in order to salvage distressed timber and alter management plans 
to minimize the impact of tree mortality on our valuable forest resources. 

Remote sensing research conducted over the past six years has attempted 
to previsual ly discriminate root-rot ( Poria wei r i i ) infected Douglas-fir 
trees from healthy trees in the Pacific Northwest. Unfortunately, the 
foliar color characteristics of healthy and diseased trees are so similar 
that aerial photographic techniques in the visible and near infrared 
portions of the spectrum (.3^ micrometer to .90 micrometer) are unable 
to differentiate affected trees. Thermal infrared detection of incipient 
root-rot infection centers was not successful with the University of Michigan's 
mul tispectral scanner Q) , although some earlier success was reported with 
a nonimaging radiometer from a low-flying heiicopter (J_ and ^) . 

Results of two years of biophysical studies at the Wind River test 
site revealed that when infrequent tree crown temperature differences occur 
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in Douglas-fir they are not necessarily due to differences in tree vigor. 
Consequently, it was not possible to separate the diseased from healthy 
trees on the basis of infrared temperature. 

To better relate ground data to airborne remote sensing, we installed 
a comprehensive network of biophysical and tree physiological instruments. 
Sensors were placed in the trees as well as over the trees (on a double 
aerial tramway system) to record data on all parameters of stress in 
healthy and diseased trees. Details of the construction and operation 
of the aerial tramway system, the collection of data through the Vidar 
5^03 digital data logger at the test site, and the analysis of the ground 
data collated with the mul t ispectral data collected by the University of 
Michigan airborne opt i ca 1 -mechan i ca 1 line scanner system are discussed in 
our 1969 and 1970 Annual Reports. More sensitive mu 1 1 i s pectra 1 data 
collecting, processing, and analyzing systems with higher resolution 
capabilities are needed to discriminate small biophysical differences 
between individual trees. 

A promising aerial survey technique has been discovered that may 
detect well-established root-rot disease centers in certain areas from 
orbital as well as suborbital altitudes. Openings in the forest canopy 
(not the trees per se) caused by root-rot diseased and fallen timber are 
proving to be positive signatures for surveying large forested areas to 
detect disease infection centers. This report discusses a new approach 
which is being developed and tested in the Douglas-f i r-hemlock types of 
of Oregon and Washington. 

CURRENT RESEARCH ACTIVITIES 

Openings in forest stands have distinctive characteristics that can 
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be frequently identified on aerial photographs and may be exploited as 
a new type of information by resource managers. Recent photo inter- 
pretation and ground examination of 1:15,840 scale photos prove conclusively 
that Por ia we i r i i root-rot disease has a distinctive signature which has 
been found in many forested areas of the Northwest. This signature 
appears as variable-sized circular openings (ranging from 50 to 700 
feet in diameter) in the form of infection centers which may be observed 
on small photo scales and possibly on orbital imagery. 

The major objectives of the Por ia we i ri i remote sensing study 
during the past year were (1) to determine the criteria for detecting root- 
rot centers or openings, (2) to investigate the distribution of signature 
indicators on existing photographs covering the Douglas-fir types in 
Oregon and Washington, (3) to cooperate with various forestry agencies 
in developing an intensive survey program on selected test sites in 
Oregon, (4) to determine the optimum photographic scale for detecting 
infection centers, and (5) to compare various types of mul t ispectra 1 
imagery, i.e., color, color infrared, and panchromatic photography, for 
discriminating disease infection centers or openings from healthy trees. 

DETECTING ROOT-ROT CENTERS 

The criteria for identifying openings in the forest canopy that can 
be attributable to the root-rot disease, Por ia we i r ri , vary with topography, 
vegetative types, elevational range, soil and site characteristics, and 
environmental factors. The interrelationships of these different charac- 
teristics to identify broad differences between openings occuring in various 
parts of the Douglas-fir subregion can only be inferred at this time. 

Further remote sensing research is contemplated to determine meaningful 
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characteristics of root-rot centers or openings in localized areas. 


The primary signature indicator of Pori a weirri root-rot was dis- 
covered in the high Cascade Mountains of Oregon on 1:15,840 scale photo- 
graphs, The signature consists of an unusual phenomenon of bare ground 
in a half-round circular shape surrounded by trees. A conglomerate of 
openings gives the appearance of a ringworm pattern (Figure 27 ). Determining 
the cause of the openings in the forest canopy involved a multidisciplinary 
group of scientists and foresters. Openings could be the result of insect 
or disease activity, unusual geomorphologi ca 1 features such as rock out- 
crops or soil and water deficient strips. The team ground-checked several 
openings in different age classes of the Douglas-f i r-hemlock type. The 
preliminary interpretation that the openings were caused by root-rot 
disease was verified by the multidisciplinary team and has been further 
substantiated by ground surveys with many plots this year. 

Not all openings in the forest canopy are caused by root-rot disease 
so that good interpretative judgments are required to identify Pori a 
openings in different parts of the Douglas-fir subregion. The striking 
appearance of the circular bare ground signature indicator in the high 
Cascades of Oregon was the primary criterion for identifying root-rot 
centers in other forested areas of the Pacific Northwest, In the early 
stages of development, a Por ia signature appears on the aerial photo as 
a circular opening (or "hole") in the forest stand with dead and down 
trees jackstrawed in the center and standing trees in var i ou s s tages of 
deterioration at the edges of the opening. Stand openings increase 
radially in size from the center of infection. Where brush and other 
vegetative species do not invade rapidly (as in the Waldo Lake Test Site 
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Figure 2 7. 1:31,680 scale color imagery shows characteristic pattern 

of Pori a wei ri i root-rot disease in forest stands of high Cascades in 
Oregon (Waldo Lake test site). At (1) a new center can be observed; 
at (2) an older enlarging center; at (3) a series of coalescing 
centers. Note how residual trees of less susceptible types (white 
and true fir, western red cedar, and white pine) fill in initial 
openings . 
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In the high Cascades), the circular bare ground signature is readily 
discerned. At lower elevations and on better timber sites, small trees 
(primarily hardwoods) and brush invade rapidly to obscure down trees 
and increase the problem of positively identifying root-rot disease 
centers. Both large and small pockets of bark beetle-killed trees in 
various stages of deterioration in the forest stand may also confuse 
interpretation. Considerable field experience is needed by the photo 
interpreter to delineate root-rot disease-caused openings. 

LOCATION OF PROBLEM AREAS 

As a first step in establishing test sites for further research and 
development of survey methodology in detecting Pori a wei ri i root-rot 
centers in the Douglas-fir subregion, it was considered desirable to 
know where (based on the general distribution of root-rot signature 
indicators) tree mortality was in progress. Existing photographs cover- 
ing approximately 30 million acres on 12 National Forests in Oregon and 
Washington were carefully scrutinized by two interpreters for the presence 
or absence of discrete signature indicators. Interpretations on some 
8,200 photographs at scales ranging from 1:10,000 to 1:15,840 indicated 
good disease indicators on four National Forests and likely indicators 
on e i ght others . 

From this broad overview of root-rot disease in the Northwest, 
three test sites were selected (Fi gure 28) for more detailed examination, 
photographic coverage, and analysis. Each test site encompasses 9 square 
miles (3-mile x 3-mile block) of pole and sawtimber stands. A variety of 
disease, stand, and climatic conditions, slope and elevations, and forest 
types would thus be sampled. 
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1. The Waldo Lake Test Site (Williatn N. F.) is in the high 
Cascades of Oregon, has some brushy areas, and moderate to easy terrain 
access i bi 1 i ty . 

2. The Olallie Lake Test Site (Mt. Hood N. F.) is in the 
Oregon Cascades, brushy with areas of dense reproduction, and moderate 
terrain accessibility. 

3. The Divide Lookout Test Site (Siuslaw N. F.) is In the 
Coast Range of Oregon, brushy with extensive understory of vine maple, 
rhododendron, and red alder, and has extremely rough terrain with deeply 
dissected slopes of 80 to 90%. 

INTERAGENCY COOPERATION 

The interest by forest resource managers in the new root-rot disease 
signature has stimulated a strong cooperative research and development 
program to evolve a survey method for detecting and evaluating the impact 
on the forest resources of the Northwest. The following Federal and 
State agencies are contributing manpower and/or financial help to assist 
In the collection of ground truth and in obtaining aerial photography: 
Forest Insect and Disease Control, Timber Management, R-6, U. S. Forest 
Service; Remote Sensing Work Unit, Pacific Southwest Forest and Range 
Experiment Station, U. S. Forest Service; Earth Resources Survey Program, 
National Aeronautics and Space Administration; Forest Disease Research, 
Pacific Northwest Forest and Range Experiment Station, U. S. Forest 
Service; and the Oregon State Board of Forestry. This interagency 
cooperative effort will increase the likelihood of developing a practical 
airborne survey sensing system in the shortest possible time. A well- 
developed survey system (orbital or suborbital) can assist in the 
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Root-rot Disease Test Sites 

Figure 28. Forest map of Oregon and Washington shows location of 
three 9~square-mi 1 e test sites -- two in the high Cascades and one 
in the Coast Range -- that are utilized for developing aerial photo 
graphic survey techniques from orbital or suborbital altitudes. 
Areas of Douglas-fir and hemlock forest types susceptible to the 
root rot disease are located in the shaded areas. 
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protection of one of our most valuable natural resources in the Pacific 
Northwest -- our National Forests. 


PROCEDURES 


SUBORBITAL MULTISPECTRAL PHOTOGRAPHY 


Preliminary interpretations of panchromatic photos at the 3 test 
sites (9 square miles each) indicated a wide range of terrain, vegetative 
conditions, and appearance of root-rot disease openings. This necessitated 
a preliminary ground evaluation of the different sites to determine the 
feasibility of testing minimum and maximum photographic scales. The 
vegetative complexity of the Coast Range test site clearly indicated 
the need for maximum ground detail to delineate root-rot openings. Conse- 
quently, photo scales larger than 1:15,840 were considered essential for 
the initial tests in the Divide Lookout test site of the Coast Range 
(Figure 29). Because of the more easily discerned signature indicators 
in the high Cascades of Oregon (Waldo Lake and Olallie Lake) much smaller 
scale photography (including orbital imagery) should be useful in detecting 
Pori a we i ri i disease centers. 

Subsequently, plans were made to secure aerial imagery at the 
respective test sites using the following scales of photography, camera 
systems, and film-filter combinations. 

1 . Photo Scales 


Divide L. 0. 
(Coast Range) 


Waldo Lake 
(Cascades) 


Olallie Lake 


1:4,000 (USFS) 
1:8,000 (USFS) 
1:15,840 (USFS) 


1:8,000 (USFS) 1:8,000 (USFS) 

1:15,840 (USFS) 1:15,840 (USFS) 
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Figure 29. 1:8,000 scale color photography shows two types of 

openings in forest stands of the Coast Range in Oregon. Plots 1 
and 2 are caused by Poria root rot. Dead Douglas-fir trees, 
killed in past years by the Douglas-fir bark beetle, remain 
standing in some openings (Plot 3). Large numbers of dead 
trees in a pocket indicate beetle activity rather than inci- 
dence of root-rot spread. 
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Photo Scales, continued 


Divide L . 0 . 

Waldo Lake 

0 1 a 1 1 i e La ke 

(Coast Range) 

(Cascades) 



1:30,000 (USFS) 

1:30,000 (USFS) 


1:60,000 (NASA) 

1:60,000 (NASA) 


1:125,000 (NASA) 

1:125,000 (NASA) 


1:250,000 (NASA) 

1:250,000 (NASA) 


2, Film-filter Combinations 

Each test site was photographed with Aero color negative film 
(2445) which provides either color or black-and-white prints, 
and with Color Infrared film (2443) using a Wratten #15 filter. 

3. Camera Systems 

The U. S. Forest Service used a Zeiss 8 1/4" focal-length camera 
for all film and scale combinations. NASA expects to use two 
RC-8 cameras (6" and 12" focal lengths) and four 70 mm cameras, 

1.57" (40 mm) focal length to 6" focal length (2 1/4" x 2 1/4" 
format) from an RB-57 or U-2 jet type aircraft. 

Multiscale photography was taken by the U. S. Forest Service using 
the PSW Forest and Range Experiment Station's Aero Commander and the 
Region 6 Cartographic Section's twin Beechcraft. Scheduling constraint 
for the NASA aircraft in 1971 prevented our obtaining small-scale photographic 
coverage. Plans were to analyze and compare the accuracy of interpreting 
several film-filter and photo scales taken at suborbital altitudes at the 
three test sites. This photo coverage has been rescheduled for NASA 
aircraft in July 1972. 
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COLLECTION OF GROUND TRUTH 


A plan to collect ground truth requires a knowledge of the general 
distribution of disease infection centers on each test site. Thus, 
criteria were developed to identify the appearance and characteristics 
of For ia we i r i i root-rot openings on aerial photographs using some pre- 
liminary ground inspections. Then, with these criteria, two photo 
interpreters scrutinized the currently available 1;15,840 scale pan- 
chromatic photos taken in 19&7-69 of each 9-square-mile area selected 
for the study. Openings and type areas suspected as indicative of root- 
rot centers were classified into two categories: good For ia signatures 

and possible For ia signatures. Each suspect area was circled and numbered 
on a frosted acetate overlay attached to each photo. Areas of timber that 
did not appear infected with For ia were also marked on the photos to serve 
as control checks. 

From these interpretation data, 25 or more areas in each of the above 
categories were selected for ground visitation on each test site. Areas to 
be checked were chosen largely on the basis of accessibility. The extremely 
adverse terrain on the Divide Lookout test site (Siuslaw N. F.) necessitated 
that plots be within a few chains of a road. Presence or absence of For i a 
root-rot was determined by the examination of several living, dead, and 
down trees for diagnostic symptoms and signs of the fungus, including ecto- 
trophic mycelium, setal hyphae, and laminated decay. Check plots were at 
least 1/5 acre in size and located at least 2 chains from known For i a centers. 
Two 2-man teams under the direction of a trained forest pathologist visited 
designated plots on the three test sites from July 6 to September 22, 1971. 


PHOTO INTERPRETATION 


Photo interpretation on this study consisted of four phases that 
ranged from a comprehensive overview for the large problem area (the 
Douglas-f i r-hemlock timber resource in the Pacific Northwest) to more 
intensive and detailed interpretations of specific plots within each of 
the three test sites. The first phase was mentioned earlier and in- 
volved the inspection of some 8,200 panchromatic photos covering approxi- 
mately 30 million acres to detect any type of signature indicators in 
the forest types of Oregon and Washington that might be caused by a 
root- rot forest disease. 

The second phase was the intensive photo interpretation on 1:15,840 
scale panchromatic photographs of three selected test sites in Oregon. 

The three test sites represent signature indicators ranging from simple 
to complex degrees of i nterpretabi 1 i ty . Two experienced photo inter- 
preters received field training on different types of openingsin forest 
areas of the high Cascades and Coast Range before interpreting For ia 
signatures on photos of the three test sites. From the large array of 
interpreted areas, more than 75 plots were chosen for field examination 
in each test site. Good accessibility was the primary concern for final 
plot selections. Photo interpreters from the two field crews greatly 
enhanced their ability to discriminate root-rot centers on the photos as 
the field visitations progressed during the summer. 

The third interpretation phase developed toward the end of the summer 
as new photo coverage on each test site with black-and-whi te, color, and 
color infrared films and at different photo scales became available. The 
two most experienced interpreters examined all test site areas, first on 
small scale and then progressing to large scale with a random selection 
from the available film types. 
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The fourth and final phase of the photo interpretation for the 1971 
season is now in progress. A more complete analysis of all film-filter 
and scale combinations is needed to make statistically valid and meaningful 
comparisons that can be utilized in establishing sound aerial photographic 
survey procedures. In order to achieve this objective, it was considered 
appropriate to devise a photo interpretation testing technique that would 
minimize or remove almost entirely the bias that is commonly found in 
multiple and sequential photographic interpretations of the same areas 
with different films and scales. 

The large number of plots selected and field-checked on each test 
site provided a sound basis for using a randomized block design with a 
factorial feature at each of the three test sites. Nine combination treat- 
ments are possible at each test site -- three types of photography and 
three photo scales. An even distribution of plots within each combination 
treatment for all three tests (dictated by total field plots) gave eight 
replications. Thus, from a total of 72 plots in each test site, (there 
were from 75 to 100 plots on each of the three test sites) eight were 
randomly selected for each film-scale combination. Photo interpretation 
will be made by five interpreters. Each of the 72 plots will be interpreted 
only once by each interpreter, thus removing a possible source of interpreter 
bias. Prior to making the random selection, the 72 plots were delineated 
and numbered on each film-scale combination at the three test sites -- a 
total of 1,9^4 plots. The interaction between photo type scale combinations 
(nine) and interpreters (five) is used as the error term for testing signi- 
ficance. Tests of significance will be performed for film types, photo scales 
and the photo-scale film interaction. 
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RESULTS 


The results of the photo interpretation are not complete at this 
time because color prints, panchromatic prints, and color IR trans- 
parencies were not available for interpreting all photo scales and test 
sites before the end of the field season. However, the results shown 
in Table 27are encouraging. Interpretation of each circled plot included 
in the field check indicated whether or not root-rot disease was present 
or absent, A few of the field plot openings in the Waldo Lake and Olallie 
Lake test sites were caused by another root rot. Pomes an nos us . .Character- 
istics of Pomes root rot in Douglas-fir are indistinguishable from Por ia 
we i r i i root rot on aerial photography and are included in the interpretation 
results. This distinction must be made on the ground. 

An inspection of the above photographic interpretation accuracy of 
rating stand openings shows that differences in photo scale do not mater- 
ially affect accuracy. Relatively small differences in accuracy are indicated 
between scales for each of the high Cascade test sites. It is apparent that, 
with the excellent signature indicators in the Waldo Lake test site area, a 
high degree of accuracy should be attainable in detecting Poria centers in 
this part of the high Cascades. Somewhat less accuracy was obtained in the 
Olallie Lake test site area. 

PUTURE RESEARCH ACTIVITIES 
DETERMINATION OP OPTIMUM PILM-SCALE COMBINATION 
Our major effort for the immediate future will be to complete the 
Interpretations (five) of the 216 selected Poria root-rot plots on the 
three test sites. Upon completion, the statistical analysis of the data 
will determine optimum combinations of films and photo scales that would 
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Table 27. Photo interpretation accuracy in detecting root-rot disease 
on three test sites in Oregon with several film types and scales of 
photography. 


Type of 
Photography 

Sea le 

Wa 1 do 
Lake 

Olallie 
■ La ke 

Divide 
L. 0. 

Pan 

1 

: 15,840 

93% 

80% 

64% 

Color 

1 

:31,680 

87% 

73% 


Color 

1 

: 15,840 


"V 

61% 

Color 

1 

:8,000 

89% 

76% 

62% 

Color 

1 

:4,000 



50% 

*'Not i 

i ncluded 

in the photographic 

coverage 
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be applicable to surveys in the Douglas-fir region or in particular test 
site areas in the Cascades or Coast Range of Oregon and Washington. Further 
remote sensing tests may be necessary if variable results for a particular 
test site occur. Preliminary data indicate good potentials for photo 
survey techniques in the high Cascades of Oregon using small-scale sub- 
orbital photography. 

Additional aerial photographic imagery will be secured over the 
three test sites from NASA RB-57 overflights, if possible, in July 1972. 
Imagery would simulate orbital and suborbital photo scales not currently 
available for analysis and projection into further survey methodology. 
Analysis of these data, coupled with photography currently on hand and 
in the interpretation and analysis stages, should provide corroboration 
on optimum film-scale combinations. This should indicate success-ratio 
likelihood for forest disease aerial survey techniques from orbital and 
suborbital altitudes. 

PHOTO SURVEY TEST IN HIGH CASCADES IN 1972 

The results of the preliminary interpretations i n the h i gh' Cascades 
of Oregon are sufficiently encouraging so that an attempt will be made in 
1972 to implement an aerial photographic trial survey in the forest types 
of the high Cascades in and adjacent to the Waldo Lake test site. An 
interagency cooperative survey development program (mentioned previously) 
is being planned for (1) administering this aerial survey, (2) designing 
the aerial survey methodology, (3) procuring, processing, and interpreting 
aerial photographic imagery, (4) collecting the ground truth according to 
specific sampling procedures, (5) assigning responsibilities to qualified 
personnel in all phases of the survey program, and (6) developing logistics 
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and setting time schedules to accomplish the aerial photographic survey in 
an expeditious manner. Details of the survey program will be developed 
during the winter period. 

IMPROVING PHOTO INTERPRETATION IN THE 
COAST RANGE TEST SITE 

Additional remote sensing research is being considered for improving 
photo interpretation capabilities in the Coast Range test site areas. A 
major problem in interpreting imagery in the Coast Range of Oregon and 
Washington is the abundance of brush and hardwood tree species which when 
fully foliaged preclude seeing details in down trees and on the ground. 
Hence, there is a need to utilize photography taken in these areas prior 
to spring growth, which is before April 15- Attempts will be made from 
mid-March to early April to determine whether a higher interpretation 
success ratio can be obtained at this time. Additional photo tests may 
also be scheduled in early October when normal leaf drop of vegetation 
(hardwoods and deciduous brush species) occurs. As a part of this 
endeavor, more detailed criteria will be investigated and developed for 
classifying root disease openings in the Coast Range forest canopy. 

These criteria would be applicable to the Oregon Coast Range and the 
Washington Olympic areas. 

ESTIMATING CHANGES IN PORIA CENTERS 
Investigations will continue on comparing older photography taken 
over the past 30 years of the Ola 1 lie Lake and Waldo Lake test sites to 
determine whether or not Por ia root rot changes shape and size rapidly, 
gradually, or remains static. More survey information is being gathered 
on the distribution of root-rot disease centers to establish a base map 
from which comparative measurements can be obtained. Disease enlargement 
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comparisons will be made on sequential photography using a Zeiss 
stereotape-Ol ivetti computer program. A time sequence for rephotography 
may also be ascertained from this investigation. 

ERTS IMAGERY MAY BE IMPACT SURVEY TOOL 
Because of the relatively large size of the root-rot infection 
centers, there is a good possibility that ERTS imagery might detect the 
incidence of advanced stages of Por ia . The three Pori a wei r i i test sites 
in Oregon are included in the Oregon State University ERTS proposal. Plans 
have been made to look at the ERTS imagery of these test sites when it 
becomes available and to collate the data with suborbital imagery taken 
in 1971 . A multistage probability sampling design may be an effective 
means to determine the Impact of the root-rot disease on the forest 
resource using ERTS imagery as the first stage. 
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CALIBRATION OF COLOR AERIAL PHOTOGRAPHY 


by 

R. W. Dana and R, J. Myhre 
INTRODUCTION 

The guiding premises behind film calibration work by this group 
are that film response data would aid in the design of remote sensing 
experiments and provide the photographer with film speed data of new 
film emulsions. The capability to relate film density to irradiance 
level and reflectivity of ground targets is also of great interest. 
Furthermore, it is reasonable to assume that calibration of this sort 
can be a valuable tool in the automatic recognition of forest and range 
type classes and their associated patterns. 

The acquisition of these capabilities requires good instrumentation -- 
a sensitive, versatile densitometer, which we have purchased, and a repeat- 
able uniform light source (sens i tome ter) which we are developing. The 
subject of a large portion of this report is the design of filter combina- 
,tions for the sensitometer that will yield a close fit to typical daylight 
spectral distributions in the visible and near infrared region, A detailed 
radiometric measurement of the sensitometer exposure time and repeatability 
was made, and a brief description of the techniques and results should be 
of interest. 

During a lull in the summer flying schedule, we were able to field- 
test some new film emulsions, using sensitometry to monitor the exposure 
and processing. Two of the film types were new to us (Eastman 24A-3 and 
SO-397) and the others were new emulsions of film types with which we had 
much previous experience. Most of the qualitative analysis can be presented. 
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while the quantitative phases are still under way, only the methods of 
analysis can be described in this report. 

SENSITOMETER DEVELOPMENT 

The instrumentation work involved with sensitometry is a continua- 
tion of work described in a previous special report (2^). The important 
requirements of an aerial film sensitometer which we have strived for 
are the fol lowing: 

1. Daylight spectral balance in the range 400-900 nanometers. 

2. Radiometric repeatability approaching the densitometer accuracy 
of + 0.02 density units. 

3. Exposure time of under 20 milliseconds. 

4. Compact size. 

5. Filter interchangeability. 

This section discusses the quantitative work toward accomplishing require- 
ments 1 , 2 and 3 . 

SPECTRAL DISTRIBUTION 

If one were to choose one particular spectral distribution for his 
light source in aerial film calibration, it should be the daylight that 
illuminates his targets. But daylight takes so many forms as evidenced 
by the variations in irradiance with wavelength that have been discovered 
by several authors. Gates {J. snd 4) has published numerous daylight spectra 
of direct sunlight from various solar zenith angles and with various atmos- 
pheric effects present. He has also presented the distribution of skylight, 
which is largely Rayleigh scattering by air molecules, and "cloudlight" due 
to overcast conditions. The latter two are much stronger in the blue end of 
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the spectrum than direct sunlight. The effect of sun angle is shown 
In Figure 30. Each curve is normalized to the same value at 480 nm to 
allow comparison of general shape. A low sun situation (large zenith 
angle and large air mass) is richer in longer wavelength radiation than 
a noon sun situation except in the very strong absorption bands. 

The sunlight spectrum is characterized by approximately 5900^K 
blackbody radiation undergoing strong absorption by ozone and oxygen 
In the ultraviolet, a moderate molecular oxygen absorption band at 
760 nanometers, and a strong water vapor band centered at 930 nm (Figure 30), 
Gates (2) has computed these attenuation effects as well as Rayleigh 
scattering by air molecules and Mie scattering by aerosols. The aerosol 
levels expected in aerial photography are rather hard to predict, and the 
effect of water vapor concentration is generally no greater than the 
variation of Figure 30. The moderate slant path represented by m = 1 .5 
would seem to be the best bet for consideration as a typical direct sunlight 
curve. 

The effect of adding skylight to the sunlight distribution is shown 
In Figure}!, revealing a general increase with decreasing wavelength in 
the 400-900 nm band. Here, as with all spectra in this section, distri- 
butions are shown as log plots in terms of relative energy per unit area- 
time-wavelength . 

Data from two other sources (6 and \ 2 ) are tabulated in Tables 28 
and 29 at 20 nm intervals along with the Gates data (3., 4). For the visible 
range covered by panchromatic and normal color films, the daylight values 
were lowered at the short wavelength end by the attenuation of a typical 
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RELATIVE SOLAR IRRADIANCE 



1+00 500 600 700 ■ 800 900 1000 

WAVELENGTH IN NANOMETERS 


Figure 30. Relative spectral irradiance at sea level from direct sun- 
light with air mass m of 1.0, 1.5, and 2.0. The corresponding sun angles 
are 0 , 48 11', and 60 from the vertical. Precipi table water concentra- 
tion -- 1 centimeter; aerosol concentration -- 200 particles per cubic 
centimeter; ozone concentration -- 0.35 centimeter. Data derived from 
reference no. 3. 
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LOG RELATIVE IRPADIANCE 


0.6 + 


O.U I 


0.2 t 


0 I 


WAVELENGTH IN NANOMETERS 


Figure 3 !• Relative spectral irradiance at sea level of direct sunlight 
at m « 1.5 and sunlight plus skylight. Data derived from references no. 3 
and no. 4. 



Table 28. Daylight spectral distributions for the visible region normalized 
to the same integrated 400 to 700 nm value. 


Wavelength 

(nm) 

DIA 

DIB 

Log Re 
D 2 

ilative Ii 
D 3 A 

'radiance 

D 3 B 

d 4 a 

D 4 b 

400 

1.53 

1.43 

1.44 

1.43 

1.67 

1.45 

1.54 

1+20 

1.64 

1.58 

1.54 

1.53 

1.68 

1.61 

1.70 

1+1+0 

1.64 

1.59 

1.60 

1.62 

1.72 

1.67 

1.71 

1+60 

1 .T 4 

’ 1.81 

1.69 

1.70 

1.76 

1.74 

1 . 75 - 

480 

1.75 

1.93 

1.72 

1.72 

1.74 

1.74 

1.75 

500 

1.72 

1.95 

1.70 

1.71 

1.72 

1.72 

1.73 

520 

1.69 

1.96 

1.68 

1.71 

•1.70 

1.70 

1.70 

540 

1.71 

1.91 

1.71 

1.72 

1.70 

1.70 

1.69 

560 

1.70 

1.83 

1.71 

1.71 

1.68 ■ 

1.71 

1.67 

580 

1.69 

1.73 

1.70 

1.71 

1.66 

1 . 68 

1.65 

600 

1.68 

1.60 

1.70 

1.69 

1.63 

1.66 

1.63 

620 

1.66 

1.54 

1.70 

1.68 

1.62 

1.65 

1.61 

640 

1.65 

1.46 

1.69 

1.68 

1.60 

1.63 

1.63 

660 

1.63 

1.39 

1.68 

1.67 

1.58 

1.62 

1.62 

680 

1. 58 

1.39 

1.63 

1.66 

1.57 

1.62 

1.62 

700 ’ 

1.60 

1.37 

1.65 

1.63 

1 

1.53 

1.60 

■rai 

Total 

25.06 

25.07 

24.99 

25.04 

24.96 

24.97 

25.05 

j 


DIA 

DIB 

D2 

D3A 

D3B 

D4A 

D4b 


sunlight + s ky 1 i ght 

cl oudl i ght 

sunl i ght m = 1 .5 

5500°K correlated color temperature 
6500°K correlated color temperature 
6000°K sunlight + skylight 
6500°K overcast skylight 
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Table 29. Daylight spectral distributions for the visibler 
infrared region normalized to the same integrated 500-920 nm 
value. 



Lc 

DIA 

)g Relati\ 
DIB 

re Irradia 
D 2 

i,nce 

D 3 A 

500 

0.83 

1.26 

0.77 

0.79 

520 

2.36 

2.83 

2.31 

2.35 

5H0 

2.57 

2.97 

2.53 

2.55 

560 

2.58 

2.91 ■ 

2. .55 

2.56 

580 

2.57 

2.81 

2.55 

2.56 

600 

2.57 

2.69 

2.55 

2.55 

620 

2.55 

2.63 

2.55 

2.54 

61+0 

2. 5I+ 

2.55 

2.54 

2.54 

660 

2.52 

2.1+8 

2.53 

2.53 

680 

2 . 1+7 

2.1+8 

2.48 

2.52 

TOO 

2 . 1+9 

2.1+6 

2.50 

2.49 

720 

2 . 1+7 

2.1+2. 

2.48 

■ 2 . 42 . 

7I+O 

2.1+6 

2.1+0 

2.47 

2.50 

T6o 

2.28 

2.23 

2.28 

2.29 

780 

2 . 1+1 

2.25 

2.43 

2.43 

800 

2.38 

2.18 

2.40 

2.40 

820 

2.38 

2.18 

2.40 

2.39 

81+0 

2.36 

2.09 

2.38 

2.37 

860 

2.3I+ 

2.01 

2 '. 37 

2.35 

880 

2.15 

1.95 

2.17 

2.16 

900 

2.10 

1.88 

2.11 

2.10 

920 

2.0I+ 

1.83 

2.06 

2.05 

Total 

i 

1+9.98 
1 

1+9. 9I+ 

1 

49.99 

50.02 


DIA - sunlight + skylight 

DIB - cloudl ight 

D2 - sunl i ght m = 1 .5 

D3A - 5500°K correlated color temperature 
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is D4 a which is further defined as having a 6000^K color temperature. 

DIB and D4B are overcast sky distribution. The color temperature of D3A 
is 5500^K while the somewhat more blue distribution D3B is for color 
temperature 6500^K. Distribution d4b, also said to have a 6500^K color 
temperature, is plotted with DIB and D3B in Figure 32 to show how daylight 
expressions can vary in the literature. d4b has more strength in the 
red and less strength in the blue than D3B. By comparison the DIB cloud- 
light is very strong in the green wavelengths. 

The filter matching technique used to obtain sensitometer filter 
combination to match daylight distribution was to start with some of the 
filters commonly used to convert tungsten light to daylight in the visible 
region and to reduce the 700-900 nm region. We computed their effect on the 
system and tried several other filters in the same series that showed from 
their own transmission curves some promise of bending the system curves 
in the right directions. The products of three major glass and gelatin 
filter suppliers were tried, and several more sources are to be investigated 
in the future. 

Spectral irradiance data for tungsten lamps at two different tempera- 
tures were digitized and combined with data for the sensitometer wedge, neutral 
density attenuators (for overall adjustment of light level for different speed 
films) and the various prospective conversion filters, again at 20 nm intervals 
Base 10 logarithms of these functions, with uncertainty + 0.01, were used to 
compile the product functions F(x) wavelength by wavelength. Equation 3 gives 
the general form of the spectral product functions. 
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LOG RELATIVE SPECTRAL IRRADIANCE 



1 + 00 . 500 600 . 700 

WAVELENGTH IN NANOMETERS 


Figure 32. Spectral distributions from three sources which might be 
used as typical curves of overcast conditions. DIB, D3B, and D4 b are 
from references no. 4, 6, and 12, respectively. 
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( 3 ) 


FU) = L(X) F(X) F (X) nF'(X), where 

n i c 

L(X) = lamp Irrad lance 

W(X) = sensltometer wedge transmissivity 

F (X) = neutral density attenuator transmissivity 
n 

HF* (x) = product of all conversion filter transmissivities. 
i c 

The logarithmic form is normalized by equation 4, 




F(X) = L(X) + W(X) + F^(X) 




M 


_ 25.00 for visible range 
50,00 for IR film range 


The eight best filter combinations for which reasonably reliable 
transmission data were available are described in Table 30 . The resultant 
distributions from equation h are displayed in Tables 31 and. 32. Combination 
FI is somewhat idealized in that a perfectly neutral wedge and general 
attenuator were assumed. This appears to be the assumption made for the 
sensitometry of Kodak Aerochrome Infrared Film 2443 and 3443 as shown in 
the manufacturer's specifications (9). While combinations Fl and F3 include 
correction filters which total 7 mm i n thickness and must be mounted 
rigidly and close to the lamps, combination F2 has the advantage of 
having 5 mm of glass filters and a gelat i n-based filter (Agfa 144 Mired), 
which can be mounted in a curved holder. The curved configuration helps 
to maintain uniformity at the film plane. Tests of uniformity with these 
filter combinations have not as yet been performed. 

Having normalized the distributions for daylight and sensltometer 
light sources for the bandwidths of interest, we found it a simple step to 
compute root-mean-square deviations between distributions. These are 
shown in Tables 33 and 34 and are representative of the closeness of fit 
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Table 30. Description of lamp-filter combinations for artificial daylight 
sources . 


i 

j 

Combination > 

1 

Lamp 

Color Temp ^ 

Filters 

Visible Region Infrared Film Region 

1+OO-TOO nm 500-920 nm 

FI 

2850 

Corning 5900 + 
Corning 3966 

FI (visible) + Kodak 
Wratten 12 

F2 

3100' 

i 

j 

Carbon vedge + Kodak 
neutral density + 
Agfa l^+U Mired cor- 
rection filter + 
Corning 3966 + Cor- 
ning 4602 

F2 (visible) + Kodak 
■ Wratten 12 

F3 

3100 

Carbon wedge + Kodak 
neutral density + 
Corning 5900 + Cor- 
ning 3966 

F3 (visible) + Kodak 
Wratten 12 

Fk 

3100 

Carbon wedge + Kodak 
neutral density + 
Kodak Wratten 38 + 
Kodak CC 20 blue 


F5 

3100 


Carbon wedge + Kodak 
neutral density + 
Corning 5900 + Cor- 
ning 4602 + Kodak 
Wratten 12 


U2 
















Table 31* Artificial daylight distributions for the visible 
region. Tungsten lamps in combination with various filters as 
described in Table 30. Normalization is for the whole 400-700 
nm band. 


Wavelength 

(nm) 

Log 

FI 

Relative 

F2 

Irradian 

F3 

ce 

f 4 

1+00 

“ 

1.53 

1.23 

1.50 

1.39 

1^20 

1.68 

1.52 

1.70 

1.66 

1|1|0 

I.7I+ 

1.61 

1.73 

1.78 

1^60 

1.75 

1.61+ 

1.68 

1.81 

1+80 

1.76 

1.70 

1.70 

1.85 

500 

1.78 


1.71 

i_^__ 

520 

1.73 

1.72 

1.70 

1.89 

5I+0 

1. 70 

1.72 

1.67 ____ 

1.87 

■ ■ 56O' 

l.jh 

■■ai 

1.80 

1.86 

580 

1.72 

1.82 

I.7I+ 

1.81 

600 

1.61+ 

1.79 

1.68 

1.73 

620 

1.58 

1. 80 

1.63 

1.62 

61+0 

I.5I1 

1.72 

1.58 

1.39 

. 660 

1.1+8 

.1.62 

1.51 

1.32 

680 

1.56 

1.53 

1.58 

1.27 

TOO 

1.59 

1.53 

1.59 

1.25 

Total 

2I+.96 

25.05 



2I+.96 

25.07 
1 


1^3 

























Table 32 . Artificial daylight distributions for the region of 
infrared film. Tungsten lamps in combination with various filters 
as described in Table 30 . Normalization is for the 500-920 nm band. 


Wavelengtlx 

(nm) 

Log 

FI 

Relative 

F 2 

Irradian 

F 3 

ce 

F 5 

500 

1.99 

0.73 

0.90 

1.06 

520 

2.57 

2.30 

2 . 1+5 

2.63 

3 h 0 

2.66 

2 . 1+9 

2.61 

2.71 

560 

2,72 

2.56 

2.70 

2.79 

580 

2.70 

2.61 

2.70 

2.79 

600 

2.63 

2.58 

2.65 

2.7I+ 

620 

2.57 

2.60 

2.60 

2.70 

61+0 

2.53 

2.52 

2.55 

2.59 

660 

2 . 1+7 

2 . 1+2 

2. k 8 

2.59 

680 

2.55 

2.33 

2.55 

2.66 

700 

2.58 

2.33 

2.56 

2.66 

720 

2.50 

2 . 1+7 

2 . 1+7 

2.57 

71+0 

2 . 1+2 

EiHI 

2,38 

2.1+5 

760 

2.31 

2.72 

2.29 

2.35 

780 

2.29 

2.68 

2.28 

2.29 

800 

2.25 

2.57 

2.26 

2.22 

820 

2.20 

2.1+1+ 

2.22 

2.11+ 

81+0 

2.18 

2.31 

2.20 

2.05 

860 

2 . 11 + 

2.19 

2.18 

1.98 

880 

2.10 

2.05 

2.16 

1.88 

900 

2 .Q 6 

1.88 

2.12 

1.75 

920 

2.12 

1.70 

2.19 

1.76 

Total 

1+9.91 

1+9.92 

1+9.95 

1+9.95 


]kk 


































































Table 33. Root mean square deviations between light sources 
FI, F2, F3, and f4 and the daylight curves DIA, D3A, D3B, DM , D4B, 
and DIB for the wavelength range 400-700 nm. The log values 
of Table 31 were used. 


1 

DIA j 

D3A 

D3B 

D4a 

d4b 

DIB 

1 

1 

i 

FI 

. o6 

1 

.09 

.05 

.06 

.06 

— 

F2 

.09 

.07 

.12 

.09 

.12 

— 

F3 

.06 

.08 

.05 

. 06 

.06 

— 

f4 

1 

1 

j 

— 

— 

— 

.09 
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Table 34. Root mean square deviations between light sources 
FI, F2, F3, and F5 and the daylight curves DIA, D2, D3A, and 
DIB for the wavelength range 500-920 nm. The log values of 
Table 32were used. 



DIA 

D2 

D3A 

DIB 

n 

.11 

.13 

.12 

— 

F2 

.15 

.16 

.16 

— 

1 

F3 

ON 

O 

.11 ■ 

.10 



F5 

— 

— 

— 

.12 
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between sens I tometer and daylight spectra. With one exception (F2 - 
D3A in Table 33)f the F3 distribution is the best fit in all comparisons 
with daylight curves, and F2 is the poorest. In Figure 33, F3 is plotted 
along with two daylight curves for the 400-700 nm band and in, Figure 34 for 
the 500-920 nm band. Daylight curves DIA and d 4A in the visible region 
are fairly well correlated as are D3A and D2 in the 500-920 nm band. 
Therefore, the latter members of these pairs are not shown in the figures. 
Figures 35 and 36 illustrate the matching of F4 and F5 with the rather hard- 
to-fit DIB cloudlight distribution. 

EXPOSURE TIME MEASUREMENTS 

The exposure time of a sens i tometer used with aerial film should be 
short enough to avoid reciprocity differences between data derived from 
sensitometry and actual imagery taken at 1/200 to 1/2000 second. To 
insure that enough light will be available for slow films and changes in 
filtration, the effective exposure time cannot be quite that fast in our 
instrument. Indications are that an adequate trade-off exists at 1/50 
second. 

The primary problem of measurement of effective exposure time is 
that the total or integrated exposure of the light pulse must be determined 
because the pulse is not simply a square step function. Total exposure E 
is defined as: 


(5) E = 7H(t) dt where H(t) is the irradiance 

o 

of the time-varying pulse. The effective exposure time T is given by 


equation 6 in terms of the maximum value H of irradiance in the pulse. 

max 
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Figure 33 . Sens 
wavelength regioi 




LOG RELATIVE SPECTRAL IRRADIANCE 



Figure 35. Cloudllght distribution, DIB, and a matching sensitometer 
spectral distribution, F4, for the visible region. 
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Figure 36. C 
spectral dist 



( 6 ) 


T _ /H(t) dt _ E 

^ H 

max max 


Thus, T is the equivalent exposure time if the pulse were a square wave 

of irradiance H . The variation of a series of E values would exhibit 
max 

the repeatability of the exposure. \ 

To integrate the pulse we imaged a 1 x 5 mm portion of the sensi- 
tometer platen onto the Photometric Data Systems microdensitometer and 
measured the relative irradiance with a 14 x 700 pm rectangular scanning 
slit. The linear amplifier used for. transmission measurements was 
adjusted to give a 0-1.2 volt output. This analog output was displayed 
on an oscilloscope and fed into the digital data system. Figure 37 is a 
scope camera photo of the pulse showing the pulse half-width to be 
approximately 20 milliseconds. 

The gearboxes on the microdensitometer film transport were disengaged 
so that a scan run could be made without physical movement of the system. 
The film transport provided a flat two-dimensional bench for the auxiliary 
optics. A scan run of about 800 milliseconds duration was performed 
coincidently with the sensitometer exposure cycle. The measuring incre- 
ment was set at the minimum value of 2 milliseconds. 

From the digital magnetic tape the^ programmer was able to pick off 
the light pulses by means of a discrimination program that separated its 
leading and trailing edges from sensor system noise. The relative irradi- 
ance was measured and recorded with an uncertainty of about H /250. 

max 

The mean value for eleven measurements of effective exposure time 
was 19.42 milliseconds. The standard deviation of integrated exposure 
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Figure 37- Oscilloscope photo of sensitometer light pulse as measured 
by the microdensitometer. The ordinate represents relative irradiance. 
The time base is 10 milliseconds for each large division. 
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divided by the mean E value of all scans was 0.038. The conversion 


of exposure variation to approximate film density variation is made possible 
by equation 7: 

(7) 9D = Y S(LogioE) , 



2.3 E 


The range of ganmia values encountered in our microdensity measure- 
ments of color reversal film is approximately 1.3 to 2.3. These values 
yield density repeatability values of 0.021 to 0.038 and are very close 
to the goal of 0.02. Part of the variation in exposure could actually 
stem from the coarseness of the time base in the pulse integration. 

Ideally, one might expect a more accurate measure if he could sum up 
the relative irradiance of a 50 ms wide pulse at 1.0 or 0.5 ms intervals. 

The design work and tests discussed in this section have brought 
the instrumental development of film calibration nearly to completion. 

At least one more infrared absorbing filter is contemplated for trial 
in matching spectral balance to daylight. More information about the 
relevance of various daylight distributions will also be sought. Finally, 
a brief photographic test of uniformity and repeatability of the sensi- 
tometer used in conjunction with in-house processing is planned. 

APPLICATION OF SENS ITOMETRIC CALIBRATION 
Is film sensitometry worth all the bother? Can aerial film calibra- 
tion be a useful aid to a remote sensing research group? These questions 
are answered affirmatively in this section as two applications of sensi- 
tometric calibrations are discussed. One deals with the testing of new film 
emulsions -- the other with the prospects of enhancing imagery for automatic 
discriminate analysis. 
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AERIAL COLOR FILM TEST 

As new aerial films emerge on the market, a few questions immediately 
come to mind. How does this new film compare with the existing film we 
have been using? How do we properly expose this new film from an aerial 
platform, and what is the best film-filter combination to use? Can we 
process this new emulsion with our existing equipment, and what are the 
processing times and procedures? 

In an attempt to answer these questions we decided to run an aerial 
film test over an area that would simulate forestry subjects. A site was 
selected near San Pablo Reservoir, about 6 miles east of Berkeley, California. 
The test area contained stands of Monterey Pine and open areas of range 
vegetation (Figure 38). 

Color panels and resolution targets were placed on the ground in 

the open range area. The targets were four-foot square panels painted 

10 

to match standard Munsell system colors. Standard U. S. Air Force bar 
targets (black bars) were painted on white panels. 

This test was not intended to be a final test but rather a preliminary 
look to gain working knowledge of the new films. A more detailed aerial 
test of new films and existing films will be run at a later date. 

On August 23, 1971, the test was flown with the Remote Sensing Work 
Unit's plane, an Aero Commander 500B. Two 70 mm Maurer KB-8 cameras with 
6- i nch focal length lenses were used. Flights were made at 800 feet above 
mean elevation to give a scale of 1:1,600. This scale was selected because 
it is often used in insect and disease evaluation. This large scale has 


Munsell Color Company. Munsell Book of Color . Ed. 1920-60. 
Baltimore, Md. 
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test area at San Pablo Reservoir, California. Print made from infrared color 












the further advantage of accommodating small, easily produced and 
easily handled test panels. 

This first test included two new emulsions of films that are presently 
in use on remote sensing studies and two films totally new to us. Speci- 
fically, they were as follows: 

Anscochrome D/200 - Film presently in use 

Ektachrome Infrared Aero (8443) - Film used in past 

Aerochrome Infrared (2443) - Replaced Type 8443 

Ektachrome EF Aerographic (SO-397) - Possible future use 

When one works with new films or new emulsion numbers, the question 
arises, "What is the speed of this film?" From past experience we have 
found that the speed can vary from one emulsion number to another and 
is not adequately specified by the manufacturer. To each of their 
aerial color and infrared-sensitive films Eastman Kodak assigns an Effec- 
tive Aerial Film Speed (EAFS). The effective speeds are determined 
empirically, usually by actual flight tests with processing under auto- 
matic, mechanized conditions. Eastman Kodak recommends determining a 
new EAFS for processing methods other than with automatic equipment. 

The EAFS is designed to be used with the Kodak Aerial Exposure Computer (8). 

For low-altitude large- to medium-scale photography, we can get the 
best exposure by reading the subject of interest with a one degree field- 
of-view light meter. But in order to use a meter, we must have an ASA 
rating that can be used on the meter calculator. 

As a starting point for determining exposure and camera settings 
for the aerial test, an arbitrary ASA rating of 200 was assigned to all 
the films in the test. On a preliminary flight over the test area a meter 
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reading was taken of dense range vegetation. Readings were also recorded 
for the average timber and average range vegetation. A normal exposure 
setting was determined based on ASA 200 and the light reading of 
dense range vegetation. To bracket the normal exposure and assure 
a proper exposure for both timber and range conditions, two additional 
camera settings were used -- one stop over the normal and one stop under 
norma 1 . 

Wratten No. 12 filters were used with the color infrared films; no 
filters were used with the color films. An additional flight was made 
with the Aerochrome IR (2443) with a Wratten 12 filter plus a Color 
Compensating Filter CC30B. The purpose of this test was to see the 
effects of altering the color balance. 

All of the films exposed on the aerial test were processed in the 
Remote Sensing Work Unit's photo lab in Berkeley, California. Color film 
processing is handled in Nikor 70 film processing equipment (Figure 39), 
which utilizes 3 1/2-gallon capacity s ta i nless s tee 1 tanks and spiral 
reels that accommodate up to 100 feet of film. Shown below are the 
chemical processes used, chemical volumes available from manufacturer, 
and the temperature recommended for each process. 


Film 

Process 

Unit Size 

Process Temperature 

Anscochrome D/200 

AR-1 

3 1/2 gal. 

68° or 75° F. 

Ektachrome IR - 8443 

E-3 

31/2 gal. 

75° F. 

Ektachrome EF - SO-397 

E-4 

3 1/2 gal. 

85° F. 

Aerochrome IR - 2443 

E-4 

3 1/2 gal. 

85° F. 


Prior to processing, sensitometer exposures were made on each film 
using a Joyce-Gevaert Type 2L sensitometer with a continuous density 
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Figure 39- Nikor 70 mm film processing equipment in the remote sensing 
photo lab. 
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wedge. The Anscochrome and Ektachrome Infrared (8443) films were 
then processed according to the manufacturers* recommended processing 
procedures . 

Aerochrome Infrared (2443) vvas processed in standard Ektachrome film 
chemicals. Process E-4. The processing of 2443 in Nikor equipment differs 
from the processing instructions that are packaged with the E-4 chemicals. 
Agitation procedures and processing times are different. The pH of the 
color developer was adjusted by adding 33 ml of 7N sodium hydroxide 
solution to the 3 1/2 gallons of developer solution. Processing times 
and agitation procedures are given in an Eastman Kodak Data Release (2.) • 

The processing of Ektachrome EF (SO-397) was not as straightforward 
as were the other film processes. Processing procedures were requested 
from Eastman Kodak in Rochester, New York, and they replied that they do 
not actually have processing instructions for Ektachrome EF (SO-397) or 
processing cycles for 70 mm Nikor- type equipment. They did suggest that 
the instructions for Ektachrome MS Aerograph ic (2448) would be a good 
starting point, but the instructions for film 2448 are recommended pro- 
cessing procedures and times for rewind processing equipment. In order 
to get some idea of what the processing times should be for SO-397> a 
list was tabulated to show the processing times by step for 2448 (rewind 
processing), 2443 (rewind processing), and 2443 (spiral reel processing). 
With the exception of the first developer, the rewind times for 2448 and 
2443 were the same for all steps. Rewind times were then compared with 
the spiral reel processing times for 2443. This comparison showed a 
time reduction in about half of the processing steps. On the basis of 
these comparisons, we felt safe in using the processing times and 
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agitation procedures for 2443 spiral reel processing {]_) . The only 
change we made was in the first developer time -- 5 1/2 minutes instead 
of 7 minutes. The 5 1/2-minute first developer time was selected by 
noting what relationships existed between times for the different 
processes, plus the fact that Nikor processing usually requires less 
development time than rewind processing. A short length of film, 
containing a series of test exposures and a continuous wedge sensito- 
meter exposure, was processed according to our selected times. The 
processed sensitometer exposure appeared to have a good range of densi- 
ties, and the test exposures had good color balance and exposure. The 
film from the aerial test was then processed in the same batch of 
chemi cals. 

An additional processing test was run to check the effects of dif- 
ferent first developer times on the SO-397. Three 10-inch lengths of 
film were cut from film bearing the same emulsion number as the film 
used on the aerial test. Identical sensitometer exposures (with a 
continuous wedge) were made on each of the three strips of film. The 
three strips were clipped to 8- x 10-inch film hangers, and the hangers 
were coded with tape for easy identification in the dark. Then the test 
strips were processed in the Nikor tanks containing E-4 processing solu- 
tions. Processing times and agitation procedures were the same as those 
used on the test film. After the first two steps (Prehardener and 
Neutralizer) the hangers were transferred to the First Developer tank. 

The three developing times selected were 4 1/2, 5 1/2, and 6 1/2 minutes. 
At the end of the required time, a hanger was removed according to its 
coded number and placed in the First Stop Bath tank. For the remainder of 
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the processing, all three strips received identical treatment. 

As one might expect, there was little change in color balance due 
to variation in first developer time but considerable effect on film 
speed. The characteristic curves using a clear scanning aperture in 
the Photometric Data Systems microdensitometer are shown for the three 
developing times in Figure 40. The speed of f i Im deve loped for 4 1/2 
minutes is roughly one half of that for 6 1/2 minutes developing time. 

Five sets of scan runs with 100 density measurements each were 
performed on the film strips for the determination of film granularity. 

The mean density of the five runs was 1.0. The spot size was 44 pm, and 

the sampling interval was 48 pm x 480 pm x-y matrix on the film. No 

D. C. drifts in the density of the x-y arrays were seen, thereby verifying 
each scan run as having a sufficiently uniform exposure for granularity 
measurements. The standard deviation of the density multiplied by 1,000 
is the usual expression of granularity. The results displayed below shown 
a tendency for increased granularity with developing time, but the dlffer- 
ences are not significant. 

Developing time (m) Mean granularity of 

SO-397 5 scan runs Standard deviation 

6.5 14.9 1.1 

5.5 15.5 1.8 

4.5 13.4 0.3 

The indication is that for greater film speed without appreciable loss 

of acuity even longer developing times might be possible. Tests of reso- 
lution as well as granularity should be performed on the longer times, 
however . 
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Figure AO, Characteristic curves for SO-397 processed at three different 
first developer times. Densitometry done with clear filter (white light) 
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After processing, all rolls of film were cut into shorter lengths 
according to exposure (normal, over, and under) or different film-filter 
combinations. A master list of all these film strips was compiled, and 
each strip was given a random number in order to eliminate as much bias 
as possible in the evaluation. The evaluation consisted of a group of 
subjective questions given to five experienced photo interpreters. The 
interpreters had a' total of 70 years experience (ranging from 3 to 25 
years) in looking at color and color infrared films on various remote 
sens i ng studies , 

The questionnaire asked each interpreter to look at various numbered 
film strips, make certain comparisons, and give an evaluation. The flow 
chart (Figure Al) illustrates the various combinations and comparisons 
included in the interpretation. The interpreters were also asked to 
evaluate the color panels on the D/200 and SO-397 for color fidelity. 

A set of miniature color panels having the same paint as the panels 
displayed on the ground was used to make the comparisons. 

An ASA rating was achieved for each film by evaluating the inter- 
preters* selections as to which film strips had the best exposures. For 
example, the interpreters were told to look at film strips #14, and 
#16 and select the best exposure for timber. In this particular case, 
three interpreters selected #4 (normal exposure) and two interpreters 
selected an exposure between #4 and #16 (between normal and over). From 
our aerial photography records we knew the camera settings and the light 
meter readings of timber and range vegetation. By putting these values 
into the meter calculator we worked backwards to a resulting ASA rating 
for a particular film-subject combination. This was done for all film- 
subject combinations looked at by each interpreter. The average ASA 


COLOR FILIC: •: 
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Figure 4l. Flow chart showing the various combinations of film and filters 
on the aerial test and the comparisons made in the interpretation analysis. 


164 




















ratings we arrived at are shown i n the following list of ASA exposure 
i nd i ces , 


Film 

T i mbe r 

Range 

D/200 

325 

300 

SO-397 

325 

280 

2443 w/#12 

350 

- 320 


Each of these ASA ratings is actually a film speed for a particular emul- 
sion number that is exposed at low altitude (800 feet) and is processed 
by the spiral reel (Nikor) equipment and procedures. > . 

The interpreters compared the color panels on the D-200 and SO-397 
films to the miniature color panels and rated the films. For overall 
color rendition and range of color sensitivity, the interpreters selected 
the SO-397- The following listing shows the composite evaluation of five 
interpreters for Munsell color rendition on the two color films. 


COLOR 

Munsell notation 

D/200 

FILM 

S0-3i 

10 R 

Best 

— 

5 YR 

Same 

Same 

10 YR 

— 

Bes t 

5 Y 

Same 

Same 

5 GY 

— 

Bes t 

10 GY 

— 

Bes t 

5 G 


Best 


Table 35 summarizes the interpreters' subjective comments on film 
response to natural subjects. They were asked to evaluate one film 
against another in terms of overall’ color balance and the color rendition 
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of timber and range based on their past experience with color and 
color infrared films. 

The characteristic curves of the two normal color films are illus- 
trated in Figures 42 and 43- They are three-color microdens i tometer scans 
of sens] tometer exposures that were included in the test film processing 
runs. The similar location of the curves along the exposure axis verifies 
that the films are approximately equivalent in speed. The slope of the 
D/200 curve is steeper in the medium density range, indicating it has 
greater sens j t i V i ty but less latitude than has SO-397* Although the D/200 
color curves are grouped very closely together, a slight yellowish red 
tendency is suggested. The SO-397 color balance is definitely cast 
towards cyan in this example, as was also seen to be true in two other 
processing runs. These color balance conclusions are somewhat verified 
by the interpreters* opinions of color panel renditions as shown in the 
preceding table. 

A scheme for quantitative color fidelity measurements has been devised 
using reflectivity measurements. Each of the color panels was subjected 
to diffuse reflection density measurements using a daylight source and 
six narrow-band analyzing filters. The same filters can be used in the 
microdensitometer, and the detectors in the two instruments are nearly 
identical in spectral response. The microdensitometer work on the film 
test images has not been completed, but it is hoped that these data can 
be correlated with the reflectivity data to rate the color accuracy of 
films. For example, a high correlation between film density and reflec- 
tion density of a given panel should indicate that a film is very capable 
of reproducing that particular hue and chroma. A low correlation between 
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Table 35* Summary of five interpreters' evaluation 
of color and color infrared films. 


TASK 

RESULTS 

GENERAL COMMENTS OF INTERPRETERS 

Best color f i Im 
for timber 

3 out of 5 said 
SO-397 (vs D/200) 

Better discrimination between shades 
of green. More color variation 
between trees . 

Best color film 
for range 
vegetat ion 

3 out of 5 said 

SO -397 (vs D/200) 

. 

Good separation between green and 
brown vegetation. Better for 
separating small range plants by 
minor color differences. 

Best color IR 
film (f i Im- 
f i 1 ter comb i n- 
ation) for 
timber and 
range 

vegetation 

A1 1 5 agreed that 
24^3 w/#l 2 f i 1 ter 
was best for over- 
all vegetation 

1 

24^3 w/#12 - Good IR response. Best 
discrimination of vegetation, 
ground cover, and soils. Good range 
of hues - good differentiation on 
trees and between plant species. 

2443 w/#12 + 3GB - Oversaturated 
with red (especially for timber). 
Strong red masks subtle hue dif- 
ferences. No blues or greens. 

Areas on 2443 w/#12 that were 
blue are now broken into reds and 
purples. 

8443 w/#12 - Weak red response, 
mostly magenta, blue, and purples. 
Fair fbr timber interpretation. 

Good hue differentiation on timber. 
IR response poorer than on 2443 
w/# 1 2 . 

Additional test 
of f i 1 ters wi th 
2 AA 3 over range 
area in Color- 
ado one week 
after film 
test in Calif- 
ornia. #12 
f i 1 ter vs #12 
+ 3 OB. 

i 

A1 1 interpreters 
felt that the 
addition of the 
bl ue f i 1 ter was 
the best for range 
cond i t ions 

2443 w/#12 - IR response poorer for 
vegetation. Better discrimination 
between living and dead vegetation. 
IR response of plants is somewhat 
masked by an overall blue. Con- 
trast between red vegetation and 
blue background not as good. 

2443 w/#12 + 3 OB - Best discrimina- 
tion by IR response. Best range of 
colors for IR response of plants 
(red, pink, yellow, brown, green). 
Better for relating plant condition 
(maturity and IR reflectance). 
Better contrast between red and the 
green background. Flowers show up 
■ better. 
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Figure h2. Color response curves of the Ansco D/200 film test emulsion 
Kodak Wrat ten filters 92,93, and Sk were used In the densitometry for 
red, green and blue response, respectively. 


MICRODENSITY 




film image densities of two different panels might indicate an ability 
to separate those Munsell colors. Of course, exposure and sensitometry 
differences between films occur in the test. Therefore, one should use 
equivalent log E rather than film density in the regression analysis. 

We are also proceeding with film and system resolution tests. The 
subject is too complex to discuss fully in this report, but it is suf- 
ficient to say that we intend to generate modulation transfer functions 
from microdensitometer edge analysis ( j_3) . Combining these functions with 
threshold modulation figures will give the most meaningful quantitative 
measure of resolution. A test of color fidelity and resolution for high- 
altitude applications is also planned for the near future, 

IMAGE ENHANCEMENT IN PATTERN RECOGNITON 
One of the. most promising aspects of the use of photographic calibra- 
tions is the possibility of enhancing photographic imagery to aid in the 
recognition patterns and in the automatic identification of forest and 
range objects. Specifically, the quantitative adjustment for exposure 
variations and for different color balance situations is often needed. 

The color balance aspect is especially important in the case of studies 
using infrared color film which has exhibited unpredictable variation in 
response. During the time period of an extended study (1) different 
processing techniques on film emulsions might be used, (2) different 
filters might be experimented with, and (3) aging changes or even manu- 
facturing changes might take place. 

The problem of exposure variations can be exemplified by considering 
the phenomena commonly called vignetting which is the radial fall-off 
of illumination on the image plane of a camera. Two effects are present 
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here. One is true vignetting — the limitation of the cone of light from 

off-axis points by the edges of the various lens elements (11). The other 

is the reduction of irradiance due to aspect angle and the increased object 

distance for off-axis points of a flat object plane imaged on a flat image 

4 

plane (^) . The latter effect' amounts to a cos 9 reduction where 0 is the 

nadir angle for vertical aerial photography. 

4 

As an example of the cos 9 effect let us take n samples of optical 
density of images of uniform targets located at each of the nadir angles 
0°, 10°, 20°, and 25°. The relative irradiance or exposure values at 
these angles are shown below; 


Nadir angle 0 (degrees) 

0 

10 

20 

25 

Cos^0 

1.00 

0.941 

0,780 

0.675 


Each group of n samples is assumed to be distributed normally about the 
same mean value D having a standard deviation of a . They could be normalized 
to the 0=0 value using sens i tometric data and lumped together as a 
group of 4n samples with a mean equal to D and the same standard devia- 
tion a . However, if these samples were combined without calibration for 
exposure difference, the distribution would not be normal as shown in 
Figure 44. The mean of this bimodal distribution is 0 + 0.15 for a 
uniform film gamma equal to 2.0 and for a equal to 0.06. The standard 
deviation is approximately 2a (or 0,12), 

The net effect of combining samples from different areas on the image 
plane is to reduce confidence in imagery that might possibly identify a 
target by mean density alone. One could conceivably have targets with a 
true standard deviation of 0.03 to 0.04 in density units and have a 
combined distribution almost as broad as shown above due to the great 
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spread in exposure values. Wherever vignetting is not precisely 
corrected by filtration there appears to be a need for accurate cali- 
bration in film pattern recognition, especially since combining measure- 
ments at different nadir angles in different sample sizes is often 
necessary and would surely result in errors for grouped mean density. 

The previous example showed the effect of exposure variations of a 
single film layer of color film or for black-and-white film. Let us now 
consider the case of infrared color film which requires consideration of 
three spectral response curves. They can be denoted as blue, green, and 
red responses which are the false color representations of green, red, 
and infrared light, respectively, incident on the film. 

In Figure ^5 we have a set of sens i tometri c curves representing the 
general case of "aim" curves for imagery acquired through the NASA Earth 
Observation Aircraft Program (10). The green curve is a thin line, and 
the tolerances of the blue and red curves with respect to the green curve 
are expressed as stippled bands. 

Let us work with a central relative log exposure level of 1.35 from 
a graybody and levels 1/2 stop underexposed and overexposed at 1.2 and 
1*5. Assume first that the blue and red curves are the upper limits (in 
density) of the tolerance bands of Figure ^5. The resultant density 
values in the three colors are shown in Table 36. The expected standard 
errors are shown in parenthesis and represent the approximate density 
reading errors one would expect. The column on the right shows the effect 
of combining equal-sized samples of data taken from photos exposed at the 
three levels 1/2 stop apart. The mean density values vary greatly from 
the high and low exposure values, and uncertainties of up to 0.39 are 
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Figure 44. Blmodal distrtbutron of density values produced by combining 
four equal-stzed normally distributed groups without cal i brat Ton . The 
displacement in mean density of the groups is due to the cos^9 reduction 
in exposure, where 0 = 0®, 10°, 20°, 25°. Film gamma is 2.0. 



DENSITY 



RELATIVE LOG EXPOSURE 


Figure 45. NASA aim curves for the general situation of Kodak 2443 Infrared 
color film. 
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indicated. Note that the mean density is not close to the density value 
for the median exposure in the case of the blue response. 

Uncertainties may be reduced considerably by normalizing the red and 
blue curves with respect to the green curve. There appears much less scat- 
ter in the data for blue-green and red-green, but it is still not down to 
the level of 0.02 that proper calibration would yield. The effects of 
underexposure seem to be more harmful to the averaging process than over- 
exposures . 

Now let us combine the effects of exposure with film response and 
processing by considering samples of different exposure and different 
sens i tometry . Assume that four measurements are taken on film having the 
same three exposure levels and having four arrangements of the blue and 
red curves with respect to the green curve. The green curve is fixed, 
and green densities are the same as i n the previous example but with 
uncertainty of zero for simplicity. The four blue-red relationships are 
high blue-high red, high blue-low red, low blue-high red, and low blue- 
low red ^21* ^22 ' matrices of Table 37), and the values 

are within the tolerances of Figure 45. Mean values of each subgroup have 
overhead lines. 

Subgroup standard deviations are large in Table 37 in terms of usual 
density measurement capabilities, and combining data for each exposure 
level again results in large standard deviations. The normalization pro- 
cess of subtracting the green density also fails to erase the uncerta i ntly . 

Although the exposure and sens i tometr i c distributions used in these 
examples are hypothetical and are not directly derived from actual experi- 
ments, the trends of optical film density are based on firm physical 
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Table 36. The effect of exposure level on the red, blue, and green dens- 
ity of infrared color film (2443). The density values are derived from 
the upper limits of the curves of Figure 44. Standard deviations are 
shown in parentheses. 



Relative Log Exposure 
1.20 1.35 1.50 

Exposure levels 
combined in equal parts 

Green density 

1. 38 

0.94 

o.6o 

0.97 


( .02) 

(.02) 

( .02) 

(.32) 

Blue density 

MEM 

0.98 

0.64 

1.05 



(.02) 

(.02) 

(.39) 

Red density 

2.30 

1.96 

1.56 

1.94 


(.02) 

(.02) 

(.02) 

(.30) 

Blue minus 

0.l6 

0.04 

0.04 

0.08 

green density 

(.03) 

(.03) 

( .03) 

(. 06 ) 

Red minus 

0.92 

1.02 

0.96 

0.97 

green density 

(.03) 

(.03) 

( . 03 ) 

(.04) 
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Table 37. The effect of exposure, film response, and processing on the 
red, blue, and green density of infrared color film (2443). The results 
of four combinations of high and low sensitometer curves for blue and red 
response with respect to a stable' green curve are shown in the 2 x 2 matri- 
ces. Standard deviations for groupings are in parentheses. Mean values 
have overhead lines. 


Green’ density 


Blue density 


Red density 


Blue minus 
green density 


Red minus 
green density 


Relative Log Exposure 
1.20 1.35 1.50 


1. 

.38 

1. 

.38 

0.94 

0.94 

0.60 

0, 

,60 

1. 

.38 

1 

.38 

0.94 

0,. 94 

0.60 

0, 

,60 

1. 

.38 

(. 

,00) 

0.94 

(.00) 

0.60 

(. 

,00) 

1. 

.54 

1, 

.54 ■ 

0.98 

0.98- 

0.64 

O', 

.64 

1. 

,12 

1, 

,12 

0.72 

0.72 

0.52 

0, 

,52 

1. 

..33 

(. 

.21) 

0. 85 

.(.13)’’ 

0.58 

(. 

,06) 

2. 

.30 

1 

.84 

1.96 

1.46 

1.56 

1, 

.02 

2. 

.30 

1, 

.84 

1.96 

1..46 

1.56 

1, 

.02 



(. 

.23) 

1.71 

(.25) 

1.29 

(. 

.27) 
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principles and could easily occur. There are certainly indications 
that film calibrations could aid in automatic interpretation studies 
such as the land-use classification program. To the extent that 
sens i tometr ic and other photometric data are available, these cali- 
brations will be applied to the Apollo 9 and RB-57 imagery. 
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